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

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

ГДИСГраницы и поздний период Тереңдетілген 21 мин оқу

Непроницаемый разлом у скважины: метод изображений

Если скважина стоит недалеко от прямолинейного непроницаемого разлома (sealing fault), давление в ней рассчитывают методом изображений (method of images): разлом заменяют воображаемой скважиной-отражением с тем же дебитом, поставленной зеркально на удвоенном расстоянии. Пока возмущение давления не дошло до разлома, скважина ведёт себя как в бесконечном пласте; потом вклад отражения растёт, и в конце концов наклон полулогарифмической прямой и уровень логарифмической производной давления удваиваются. По времени и форме этого перехода находят расстояние до разлома.

Не үшін қажет

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

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

  • Отражение с тем же дебитом на расстоянии 2L2L делает линию разлома линией нулевого перетока. Давление в скважине — сумма её собственного отклика и вклада отражения meE1(ϕμctL2/(kt))m_eE_1\big(\phi\mu c_tL^2/(kt)\big), me=qBμ/(4πkh)m_e = qB\mu/(4\pi kh).
  • Сначала скважина «не знает» о разломе: первая прямая и первая полка производной дают khkh и скин как обычно. В конце вклад отражения сам растёт как логарифм времени, и наклон прямой и уровень производной удваиваются.
  • Производная поднимается с mem_e до 2me2m_e по закону 1+e−L2/(ηt)1 + \e^{-L^2/(\eta t)}, η=k/(ϕμct)\eta = k/(\phi\mu c_t): переход виден раньше, чем в самом давлении, и тянется около 2,7 цикла (от 1 до 99 %).
  • Расстояние — по времени пересечения прямых, L=ktx/(γϕμct)L = \sqrt{kt_x/(\gamma\phi\mu c_t)}, или по любой точке ступени производной, L=ηt ln⁡[me/(Δp′−me)]L = \sqrt{\eta t\,\ln[m_e/(\Delta p' - m_e)]}.
  • Пересечение по ещё не сформировавшейся второй прямой занижает расстояние; прямая с теоретическим наклоном 2m2m через последний замер точнее.
  • Разлом в нескольких метрах от скважины скрыт влиянием ствола: видна только удвоенная полка, и без знания о разломе проницаемость получают вдвое заниженной, а скин — ложно малым или отрицательным.
Негізгі қатынастар
Δp=me[ln⁡tD+0,809+2S+E1 ⁣(LD2tD)]\Delta p = m_e\left[\ln t_D + 0{,}809 + 2S + E_1\!\left(\dfrac{L_D^2}{t_D}\right)\right]

давление в скважине у разлома после пуска; LD=L/rwL_D = L/r_w

Δp′=me[1+e−L2/(ηt)]\Delta p' = m_e\left[1 + \e^{-L^2/(\eta t)}\right]

ступень логарифмической производной; η=k/(ϕμct)\eta = k/(\phi\mu c_t)

L=0,0447k txϕμctL = 0{,}0447\sqrt{\dfrac{k\,t_x}{\phi\mu c_t}}

по времени пересечения прямых: м, мД, ч, мПа·с, 1/МПа (oilfield — 0,01217, ft)

L=0,0596k tϕμctln⁡meΔp′−meL = 0{,}0596\sqrt{\dfrac{k\,t}{\phi\mu c_t}\ln\dfrac{m_e}{\Delta p' - m_e}}

по точке на ступени производной (oilfield — 0,01624)

Lmax⁡=0,0349qB tpϕcth δpL_{\max} = 0{,}0349\sqrt{\dfrac{qB\,t_p}{\phi c_t h\,\delta p}}

дальше этого кривая восстановления давления (КВД) после tpt_p часов работы разлом не покажет при пороге δp\delta p: м³/сут, ч, 1/МПа, м, МПа

Жиі қателіктер
  • Считать любое удвоение наклона разломом: похожую картину дают и пласт, проводимость которого вдали от скважины вдвое меньше, и ложные «прямые» на участке влияния ствола.
  • Проводить вторую прямую по точкам, где производная ещё не дошла до удвоенного уровня, и получать заниженное расстояние.
  • Обрабатывать запись скважины у самого разлома как обычный радиальный приток.
  • Принимать p∗p^* первой прямой КВД у разлома за начальное пластовое давление: он занижен на вклад отражения.
  • Ожидать, что длинная КВД обязательно покажет далёкий разлом: предел задают отобранный перед остановкой объём и разрешение манометра.

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

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

Разлом как граница пласта

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

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

Почему разлом бывает непроницаемым
Механизм
Контакт с непроницаемой породой
Что происходит
напротив коллектора — глина, соль или плотная порода
Когда вероятен
смещение больше толщины пласта
Механизм
Глинистая примазка
Что происходит
глина, размазанная по сместителю, разделяет проницаемые породы
Когда вероятен
в смещённой толще много глинистых пластов
Механизм
Катаклаз
Что происходит
зёрна раздроблены, поры в зоне разлома мелкие
Когда вероятен
чистые песчаники, особенно на больших глубинах
Механизм
Цементация
Что происходит
минералы, выпавшие из растворов, запечатали поры зоны разлома
Когда вероятен
разломы, по которым когда-то двигались подземные воды

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

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

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

Скважина-отражение

Метод изображений разобран в уроке о суперпозиции: за разлом на таком же расстоянии ставят воображаемую скважину с той же историей дебита — скважину-отражение. Две одинаковые скважины симметричны относительно линии разлома, поэтому поперёк неё нет перепада давления, а значит, нет и перетока. Здесь нас интересует, как это выглядит в записи давления самой скважины.

Скважина у непроницаемого разлома: безразмерная форма
pD=12(ln⁡tD+0,809)+S⏟скважина+12 E1 ⁣(LD2tD)⏟отражение,LD=Lrwp_D = \underbrace{\tfrac12\left(\ln t_D + 0{,}809\right) + S}_{\text{скважина}} + \underbrace{\tfrac12\,E_1\!\left(\frac{L_D^2}{t_D}\right)}_{\text{отражение}},\qquad L_D = \frac{L}{r_w}
мұндағы:
  • pDp_Dбезразмерное снижение забойного давления 2πkh Δp/(qBμ)2\pi kh\,\Delta p/(qB\mu), со скином — как в уроке Решение для линейного источника и функция давления, —
  • tDt_Dбезразмерное время kt/(ϕμctrw2)kt/(\phi\mu c_tr_w^2), —
  • L, LDL,\ L_Dрасстояние от скважины до разлома и оно же в радиусах скважины, м, —
  • SSскин-фактор реальной скважины; у отражения скина нет, —
  • E1E_1интегральная показательная функция (функция E₁), —

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

Мына жағдайда орындалады
  • разлом — прямая непроницаемая граница через всю толщину однородного пласта; других границ возмущение не достигло;
  • для самой скважины — логарифмическое приближение (tD>25t_D > 25), влияние ствола закончилось;
  • отражение — линейный источник на расстоянии 2L2L; аргумент E1E_1 равен (2LD)2/(4tD)=LD2/tD(2L_D)^2/(4t_D) = L_D^2/t_D.

Аргумент вклада отражения можно записать через пьезопроводность η=k/(ϕμct)\eta = k/(\phi\mu c_t): LD2/tD=L2/(ηt)L_D^2/t_D = L^2/(\eta t). Значит, всё зависит от одного числа — времени в единицах L2/ηL^2/\eta. Пока ηt\eta t мало по сравнению с L2L^2, вклад отражения ничтожен; когда ηt\eta t много больше L2L^2, функция E1E_1 сама становится логарифмом, и pDp_D растёт как ln⁡tD\ln t_D — вдвое быстрее, чем у одиночной скважины.

Басқа бірліктер жүйелерінде (1)
Метрические промысловые единицы
Δp=0,933 qBμkh[ln⁡3,553⋅10−3 ktϕμctrw2+0,809+2S+E1 ⁣(L2ηt)],η=3,553⋅10−3 kϕμct\Delta p = 0{,}933\,\frac{qB\mu}{kh}\left[\ln\frac{3{,}553\cdot10^{-3}\,kt}{\phi\mu c_tr_w^2} + 0{,}809 + 2S + E_1\!\left(\frac{L^2}{\eta t}\right)\right],\quad \eta = \frac{3{,}553\cdot10^{-3}\,k}{\phi\mu c_t}

qq — м³/сут, μ\mu — мПа·с, kk — мД, hh, rwr_w и LL — м, ctc_t — 1/МПа, tt — ч; Δp=pi−pwf\Delta p = p_i - p_{wf} — МПа; η\eta — м²/ч.

0,933=1/(4π⋅0,0853)0{,}933 = 1/(4\pi\cdot0{,}0853); 3,553⋅10−33{,}553\cdot10^{-3} переводит время в часах в безразмерное (Анализ кривой восстановления давления).

Скважина у разлома и скважина-отражение в планеразлом1: r < L2: r = L3: r = 2Lскважинаотражениетот же дебитL2L

Схема в плане. Слева скважина, правее — вертикальная линия разлома со штриховкой, за разломом на таком же расстоянии — скважина-отражение с подписью «тот же дебит». Внизу размерные линии: L от скважины до разлома и 2L от скважины до отражения. Вокруг скважины две штриховые окружности: малая, «1: r < L», не доходит до разлома; большая, «2: r = L», касается его. Пунктирная дуга «3: r = 2L» — часть окружности вокруг отражения, проходящая через скважину.

СызбаСкважина у разлома и её отражение в плане (схема, не в масштабе). Окружности 1 и 2 — радиус исследования 4ηt\sqrt{4\eta t} вокруг скважины: пока он меньше LL, скважина работает как в бесконечном пласте; когда он доходит до разлома, производная начинает подниматься. Дуга 3 — такой же радиус вокруг отражения: когда он равен 2L2L, возмущение от отражения дошло до скважины и вклад разлома заметен уже в самом давлении; затем наклон удваивается. Линии показаны только по эту сторону разлома.
Стадии записи скважины у разлома
Стадия
Бесконечный пласт
ηt/L2\eta t/L^2
меньше 0,2
Производная
первая полка mem_e
Полулогарифмический график
первая прямая, наклон mm
Стадия
Начало перехода
ηt/L2\eta t/L^2
0,2–0,45
Производная
выше полки на 1–10 %
Полулогарифмический график
точки чуть выше первой прямой
Стадия
Переход
ηt/L2\eta t/L^2
0,45–10
Производная
от 1,1 me1{,}1\,m_e до 1,9 me1{,}9\,m_e
Полулогарифмический график
кривая, прямой нет
Стадия
Удвоенный наклон
ηt/L2\eta t/L^2
больше 10–25
Производная
2me2m_e с точностью 5 % и лучше
Полулогарифмический график
вторая прямая, наклон 2m2m

Границы стадий — по формуле ступени производной в разделе ниже: при ηt/L2\eta t/L^2 = 0,217 производная выше полки на 1 %, при 0,434 — на 10 %, при 9,5 — прошла 90 % ступени, при 25 — 96 %. Время пересечения двух прямых соответствует ηt/L2=1,781\eta t/L^2 = 1{,}781.

Две прямые на полулогарифмическом графике

Когда ηt\eta t много больше L2L^2, функцию E1E_1 можно заменить логарифмом, как для самой скважины (логарифмическое приближение). Тогда запись описывают две прямые:

Первая и вторая полулогарифмические прямые
первая:Δp=me[ln⁡tD+0,809+2S]вторая:Δp=2me[ln⁡tD+0,809+S−ln⁡2Lrw]\begin{aligned}\text{первая:}\quad &\Delta p = m_e\left[\ln t_D + 0{,}809 + 2S\right]\\ \text{вторая:}\quad &\Delta p = 2m_e\left[\ln t_D + 0{,}809 + S - \ln\frac{2L}{r_w}\right]\end{aligned}
мұндағы:
  • Δp\Delta pснижение забойного давления pi−pwfp_i - p_{wf}, Па
  • mem_eнаклон первой прямой на единицу натурального логарифма, qBμ/(4πkh)qB\mu/(4\pi kh), Па
  • tD, St_D,\ Sбезразмерное время и скин-фактор, —
  • L, rwL,\ r_wрасстояние до разлома и радиус скважины, м

Бірліктер: любая согласованная система; логарифм натуральный

Мына жағдайда орындалады
  • первая прямая — пока ηt/L2<0,2\eta t/L^2 < 0{,}2 и влияние ствола уже закончилось;
  • вторая — при ηt/L2≫1\eta t/L^2 \gg 1; наклон отличается от 2m2m меньше чем на 5 % только после ηt/L2≈10\eta t/L^2 \approx 10;
  • скин входит во вторую прямую один раз: перепад на скине есть только у реальной скважины.

Наклон второй прямой вдвое больше — 2m2m на десятичный цикл, — а её положение зависит от расстояния до разлома. Поэтому по первой прямой, как обычно, находят khkh и скин (Анализ кривой падения давления), а вторая служит проверкой и даёт расстояние. Если те же khkh и скин, подставленные во вторую прямую вместе с LL, воспроизводят её положение, картина согласована.

Приравняв две прямые, получаем момент их пересечения — время пересечения (intersection time):

Расстояние до разлома по времени пересечения прямых
tx=γ ϕμctL2k⇒L=k txγ ϕμct,γ=1,781t_x = \frac{\gamma\,\phi\mu c_tL^2}{k}\quad\Rightarrow\quad L = \sqrt{\frac{k\,t_x}{\gamma\,\phi\mu c_t}},\qquad \gamma = 1{,}781
мұндағы:
  • txt_xвремя пересечения продолжений первой и второй прямых, с
  • LLрасстояние до разлома, м
  • γ\gammaэкспонента постоянной Эйлера, —
  • k, ϕ, μ, ctk,\ \phi,\ \mu,\ c_tпроницаемость (по первой прямой), пористость, вязкость, общая сжимаемость, м², —, Па·с, 1/Па

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

Мына жағдайда орындалады
  • обе прямые надёжно выделены; вторая проведена по точкам, где наклон действительно удвоился;
  • ϕ\phi и ctc_t известны независимо: их ошибка переходит в LL как квадратный корень.

В момент пересечения ηtx/L2=γ\eta t_x/L^2 = \gamma, то есть радиус исследования 4ηt\sqrt{4\eta t} равен 2,67 L2{,}67\,L. Производная к этому моменту уже на 57 % прошла путь от первой полки ко второй: пересечение — середина перехода, а не его начало.

Басқа бірліктер жүйелерінде (2)
Метрические промысловые единицы
L=0,0447k txϕμctL = 0{,}0447\sqrt{\frac{k\,t_x}{\phi\mu c_t}}

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

0,0447=3,553⋅10−3/1,7810{,}0447 = \sqrt{3{,}553\cdot10^{-3}/1{,}781}.

Промысловые единицы (oilfield)
L=0,01217k txϕμctL = 0{,}01217\sqrt{\frac{k\,t_x}{\phi\mu c_t}}

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

0,01217=2,637⋅10−4/1,7810{,}01217 = \sqrt{2{,}637\cdot10^{-4}/1{,}781}, где 2,637⋅10−42{,}637\cdot10^{-4} — множитель безразмерного времени в этих единицах.

Разберём условную запись. Разведочную скважину пустили в пробную эксплуатацию с постоянным дебитом 110 м³/сут; глубинный манометр записывал давление около 400 ч. По сейсмике в нескольких сотнях метров от скважины предполагается разлом.

Снижение давления при пробной эксплуатации (условные данные, выборка)
tt
36 с
Δp\Delta p, МПа
0,5654
Δp′\Delta p', МПа
—
tt
6 мин
Δp\Delta p, МПа
0,7368
Δp′\Delta p', МПа
0,0537
tt
30 мин
Δp\Delta p, МПа
0,8200
Δp′\Delta p', МПа
0,0505
tt
1,00 ч
Δp\Delta p, МПа
0,8547
Δp′\Delta p', МПа
0,0501
tt
2,00 ч
Δp\Delta p, МПа
0,8892
Δp′\Delta p', МПа
0,0499
tt
5,01 ч
Δp\Delta p, МПа
0,9352
Δp′\Delta p', МПа
0,0507
tt
10,0 ч
Δp\Delta p, МПа
0,9714
Δp′\Delta p', МПа
0,0560
tt
20,0 ч
Δp\Delta p, МПа
1,0134
Δp′\Delta p', МПа
0,0668
tt
39,8 ч
Δp\Delta p, МПа
1,0636
Δp′\Delta p', МПа
0,0784
tt
100 ч
Δp\Delta p, МПа
1,1417
Δp′\Delta p', МПа
0,0895
tt
200 ч
Δp\Delta p, МПа
1,2053
Δp′\Delta p', МПа
0,0941
tt
398 ч
Δp\Delta p, МПа
1,2714
Δp′\Delta p', МПа
—

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

Полулогарифмический график пробной эксплуатации скважины у разлома с двумя прямыми0,010,111010010000,500,751,001,25Время после пуска t, чСнижение давления Δp, МПаоранжевая — первая прямаябирюзовая — хорда последних замеровсерая — наклон ровно 2m30 ч38 ч

График снижения давления от 0,5 до 1,4 МПа в зависимости от времени от 0,01 до 1000 часов на логарифмической шкале. Точки до 0,3 часа круто поднимаются, от 0,5 до 5 часов лежат на первой штриховой прямой, затем отходят от неё вверх и к 400 часам доходят до 1,27 МПа. Вторая штриховая прямая, вдвое круче первой, проходит через последние точки и пересекает первую около 30 часов; серая сплошная прямая с наклоном ровно вдвое больше первого пересекает её около 38 часов.

ЕсептеуПолулогарифмический график записи (точки). Штриховые прямые — первая (через замеры 0,5 и 5 ч) и хорда двух последних замеров из примера ниже; серая сплошная — прямая с наклоном 2m2m через последний замер. Отмечены точки их пересечения с первой прямой.
Есептеу мысалыПроницаемость, скин и расстояние до разлома по двум прямым
Бастапқы деректер (шартты)
Скважина и пласт
qq = 110 м³/сут; BB = 1,25; μ\mu = 0,65 мПа·с; hh = 14 м; ϕ\phi = 0,19; ct=1,4⋅10−3c_t = 1{,}4\cdot10^{-3} МПа⁻¹; rwr_w = 0,1 м
Первая прямая
замеры 0,501 ч (0,8200 МПа) и 5,01 ч (0,9352 МПа) — ровно цикл; производная между ними почти постоянна
Последние замеры
199,5 ч — 1,2053 МПа; 398,1 ч — 1,2714 МПа (ровно 0,3 цикла)
Шешім
  1. Наклон первой прямой: m=0,9352−0,8200=0,1152m = 0{,}9352 - 0{,}8200 = 0{,}1152 МПа на цикл. qBμ=110⋅1,25⋅0,65=89,375qB\mu = 110\cdot1{,}25\cdot0{,}65 = 89{,}375; kh=2,149⋅89,375/0,1152=1667kh = 2{,}149\cdot89{,}375/0{,}1152 = 1667 мД·м, k=119,1k = 119{,}1 мД.
  2. Первая прямая через 1 ч: Δp1ч=0,8200+0,3⋅0,1152=0,8546\Delta p_{1\text{ч}} = 0{,}8200 + 0{,}3\cdot0{,}1152 = 0{,}8546 МПа. Скин: lg⁡kϕμctrw2=7,838\lg\dfrac{k}{\phi\mu c_tr_w^2} = 7{,}838; S=1,151 (0,8546/0,1152−7,838+2,098)=1,93S = 1{,}151\,(0{,}8546/0{,}1152 - 7{,}838 + 2{,}098) = 1{,}93.
  3. Наклон по двум последним замерам: (1,2714−1,2053)/0,3=0,2203(1{,}2714 - 1{,}2053)/0{,}3 = 0{,}2203 МПа на цикл — в 1,91 раза больше первого. Удвоение почти состоялось: это признак одного непроницаемого разлома.
  4. Пересечение: первая прямая Δp=0,8546+0,1152lg⁡t\Delta p = 0{,}8546 + 0{,}1152\lg t, хорда Δp=0,6986+0,2203lg⁡t\Delta p = 0{,}6986 + 0{,}2203\lg t; lg⁡tx=(0,8546−0,6986)/(0,2203−0,1152)=1,484\lg t_x = (0{,}8546 - 0{,}6986)/(0{,}2203 - 0{,}1152) = 1{,}484, tx=30,5t_x = 30{,}5 ч.
  5. Расстояние: ϕμct=0,19⋅0,65⋅1,4⋅10−3=1,729⋅10−4\phi\mu c_t = 0{,}19\cdot0{,}65\cdot1{,}4\cdot10^{-3} = 1{,}729\cdot10^{-4}; L=0,0447119,1⋅30,51,729⋅10−4=205L = 0{,}0447\sqrt{\dfrac{119{,}1\cdot30{,}5}{1{,}729\cdot10^{-4}}} = 205 м.
  6. Вариант с теоретическим наклоном: прямая Δp=0,6724+0,2304lg⁡t\Delta p = 0{,}6724 + 0{,}2304\lg t с наклоном 2m=0,23042m = 0{,}2304 через последний замер. Разность наклонов равна mm, поэтому lg⁡tx=(0,8546−0,6724)/0,1152=1,582\lg t_x = (0{,}8546 - 0{,}6724)/0{,}1152 = 1{,}582, tx=38,2t_x = 38{,}2 ч, L=229L = 229 м.
Жауап: k≈119,1k \approx 119{,}1 мД, скин около 1,93; расстояние до разлома — 205 м по хорде последних замеров и 229 м по прямой с удвоенным наклоном.

Бұл нені білдіреді. Запись рассчитана по модели с kk = 120 мД, SS = 2 и разломом в 230 м. Первая прямая дала проницаемость и скин с погрешностью 0,8 % и 0,07 (проницаемость чуть занижена: на краях отрезка ещё сказывается ствол и уже начинает сказываться разлом), потому что проведена раньше, чем разлом заметно повлиял на давление. Расстояние по хорде занижено на 11 %: к концу записи наклон ещё не удвоился, хорда положе настоящей второй прямой и пересекает первую раньше. Прямая с теоретическим наклоном 2m2m через последний замер ошибается гораздо меньше; здесь ошибка почти нулевая, в среднем для записи такой длины она около 3 % (таблица ниже). Данные примера условные.

Ошибка расстояния по времени пересечения в зависимости от длины записи (расчёт без влияния ствола)
Конец записи, ηt/L2\eta t/L^2
5
Отношение наклонов
1,75
LL по хорде / истинное
0,74
LL по наклону 2m2m / истинное
0,91
Конец записи, ηt/L2\eta t/L^2
10
Отношение наклонов
1,87
LL по хорде / истинное
0,83
LL по наклону 2m2m / истинное
0,95
Конец записи, ηt/L2\eta t/L^2
25
Отношение наклонов
1,94
LL по хорде / истинное
0,91
LL по наклону 2m2m / истинное
0,98
Конец записи, ηt/L2\eta t/L^2
50
Отношение наклонов
1,97
LL по хорде / истинное
0,94
LL по наклону 2m2m / истинное
0,99
Конец записи, ηt/L2\eta t/L^2
100
Отношение наклонов
1,99
LL по хорде / истинное
0,97
LL по наклону 2m2m / истинное
1,00

Хорда проведена через замеры в конце записи, отстоящие вдвое по времени; прямая с наклоном 2m2m — через последний замер; первая прямая точная. Для записи урока ηt/L2≈18,6\eta t/L^2 \approx 18{,}6, и такой же расчёт даёт для неё ошибку хорды около 12 %, а прямой с наклоном 2m2m — около 3 %.

Ступень производной

На производной разлом виден раньше и яснее. Производная функции E1(a/t)E_1(a/t) по ln⁡t\ln t равна e−a/t\e^{-a/t} (производная у прямолинейной границы), поэтому после окончания влияния ствола производная описывается простой формулой, и из неё сразу выражается расстояние:

Ступень производной у непроницаемого разлома и расстояние по одной её точке
Δp′=me[1+exp⁡(−L2ηt)]⇒L=η t ln⁡meΔp′−me\Delta p' = m_e\left[1 + \exp\left(-\frac{L^2}{\eta t}\right)\right]\quad\Rightarrow\quad L = \sqrt{\eta\,t\,\ln\frac{m_e}{\Delta p' - m_e}}
мұндағы:
  • Δp′\Delta p'логарифмическая производная в момент tt на ступени, me<Δp′<2mem_e < \Delta p' < 2m_e, Па
  • mem_eуровень первой полки производной, qBμ/(4πkh)qB\mu/(4\pi kh), Па
  • η\etaпьезопроводность k/(ϕμct)k/(\phi\mu c_t), м²/с
  • L, tL,\ tрасстояние до разлома и время с момента пуска, м, с

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

Мына жағдайда орындалады
  • однородный пласт, один прямолинейный непроницаемый разлом;
  • влияние ствола закончилось; mem_e и kk известны по первой полке;
  • точки берут в середине ступени, где Δp′/me\Delta p'/m_e от 1,2 до 1,7: в начале и в конце ступени малые ошибки производной сильно меняют LL.

Ступень — одна и та же кривая для любого разлома, если время измерять в единицах L2/ηL^2/\eta. Поэтому каждая её точка даёт расстояние, и вторую прямую ждать не нужно. Середина ступени (Δp′=1,5 me\Delta p' = 1{,}5\,m_e) наступает при ηt=L2/ln⁡2\eta t = L^2/\ln 2.

Басқа бірліктер жүйелерінде (2)
Метрические промысловые единицы
L=0,0596k tϕμct ln⁡meΔp′−meL = 0{,}0596\sqrt{\frac{k\,t}{\phi\mu c_t}\,\ln\frac{m_e}{\Delta p' - m_e}}

LL — м, kk — мД, tt — ч, μ\mu — мПа·с, ctc_t — 1/МПа; Δp′\Delta p' и mem_e — в одних единицах.

0,0596=3,553⋅10−30{,}0596 = \sqrt{3{,}553\cdot10^{-3}}.

Промысловые единицы (oilfield)
L=0,01624k tϕμct ln⁡meΔp′−meL = 0{,}01624\sqrt{\frac{k\,t}{\phi\mu c_t}\,\ln\frac{m_e}{\Delta p' - m_e}}

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

0,01624=2,637⋅10−40{,}01624 = \sqrt{2{,}637\cdot10^{-4}}.

Производная записи у разлома, отнесённая к уровню первой полки, и ступень по формуле0,111010010001,01,52,0Время после пуска t, чΔp′ / mₑ10 %50 %tₓ90 %

График отношения производной к уровню первой полки от 0,8 до 2,2 в зависимости от времени от 0,1 до 1000 часов на логарифмической шкале. Точки от 0,1 до 0,5 часа лежат чуть выше 1 — это конец влияния ствола; от 1 до 4 часов они лежат на уровне 1; затем поднимаются вдоль плавной S-образной линии и к 250 часам доходят до 1,9; линия продолжает подниматься к 2. Отмечены точки на линии: 1,1 при 9 часах, 1,5 при 31 часе, 1,57 при 38 часах и 1,9 при 204 часах. Штриховые горизонтали — 1 и 2.

ЕсептеуПроизводная записи урока, делённая на уровень первой полки (точки), и ступень по формуле выше для разлома в 230 м (линия); шкала времени логарифмическая. Отмечены подъём на 10 %, середина ступени, время пересечения прямых и 90 % ступени.
Есептеу мысалыРасстояние до разлома по одной точке ступени
Бастапқы деректер (шартты)
Первая полка производной
среднее от 1 до 3,2 ч — mem_e = 0,04996 МПа
Точка на ступени
через 25,12 ч производная 0,07068 МПа
Пьезопроводность по первой прямой
η=3,553⋅10−3⋅119,1/(0,19⋅0,65⋅1,4⋅10−3)=2447\eta = 3{,}553\cdot10^{-3}\cdot119{,}1/(0{,}19\cdot0{,}65\cdot1{,}4\cdot10^{-3}) = 2447 м²/ч
Шешім
  1. Отношение к полке: Δp′/me=0,07068/0,04996=1,4147\Delta p'/m_e = 0{,}07068/0{,}04996 = 1{,}4147 — середина ступени, удобная точка.
  2. ln⁡meΔp′−me=ln⁡11,4147−1=0,880\ln\dfrac{m_e}{\Delta p' - m_e} = \ln\dfrac{1}{1{,}4147 - 1} = 0{,}880.
  3. L=ηt⋅0,880=2447⋅25,12⋅0,880=233L = \sqrt{\eta t\cdot0{,}880} = \sqrt{2447\cdot25{,}12\cdot0{,}880} = 233 м.
  4. Чувствительность: если производная прочитана выше на 2 % от полки, отношение 1,4347 и LL = 226 м — ошибка 2,7 %.
  5. Та же формула через 6,31 ч, где подъём только начался (Δp′/me=1,0366\Delta p'/m_e = 1{,}0366), даёт LL = 226 м, но ошибка производной на те же 2 % от полки дала бы уже 211 м — на 7 % меньше.
Жауап: По точке в середине ступени расстояние до разлома около 233 м; модель — 230 м.

Бұл нені білдіреді. Ступень производной содержит расстояние в каждой точке, поэтому разлом можно найти, не дожидаясь второй прямой: здесь хватило бы 40 ч записи вместо 400 (производную в точке 25 ч считают по замерам до 40 ч). Платой служит зависимость от модели: формула верна для однородного пласта с одной прямолинейной границей, а уровень первой полки нужно знать точно. На практике то же делает подбор модели с отражением ко всей записи: программа сравнивает расчётную ступень с замеренной сразу во всех точках (Эталонные кривые). Данные примера условные.

Есептеу мысалыСколько записывать, чтобы увидеть разлом
Бастапқы деректер (шартты)
Пласт
как в примерах выше; η\eta ≈ 2447 м²/ч
Разлом по сейсмике
в 230 м или в 460 м от скважины
Шешім
  1. Масштаб времени L2/ηL^2/\eta: 21,6 ч для 230 м и 86,5 ч для 460 м.
  2. Подъём производной на 10 % (ηt/L2\eta t/L^2 = 0,434): через 9,4 и 38 ч.
  3. Середина ступени (ηt/L2=1/ln⁡2=1,443\eta t/L^2 = 1/\ln 2 = 1{,}443): через 31 и 125 ч.
  4. Наклон удвоился с точностью 5 % (ηt/L2\eta t/L^2 ≈ 10): через 216 и 865 ч — 36 суток.
Жауап: Чтобы уверенно увидеть разлом в 230 м, хватит 31 ч записи (середина ступени); для второй прямой нужно около 220 ч. Для разлома в 460 м все сроки вчетверо больше.

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

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

Чем ближе разлом, тем раньше начинается переход. Если LL — всего несколько метров или десятков метров, первая полка может уложиться в минуты и целиком оказаться под влиянием ствола. Тогда в записи видна только одна полка — удвоенная — и одна прямая с наклоном 2m2m: флюид подходит к скважине только с одной стороны, как к скважине на краю полуплоскости. Такой режим называют полурадиальным притоком (hemiradial flow).

Полурадиальный приток: проводимость и скин при известном расстоянии до разлома
kh=qBμ2π Δp2′,S=Δp(t)Δp2′−ln⁡ktϕμctrw2−0,809+ln⁡2Lrwkh = \frac{qB\mu}{2\pi\,\Delta p_2'},\qquad S = \frac{\Delta p(t)}{\Delta p_2'} - \ln\frac{kt}{\phi\mu c_tr_w^2} - 0{,}809 + \ln\frac{2L}{r_w}
мұндағы:
  • Δp2′\Delta p_2'уровень удвоенной полки производной, 2me=qBμ/(2πkh)2m_e = qB\mu/(2\pi kh), Па
  • Δp(t)\Delta p(t)снижение давления в момент tt на удвоенной полке, Па
  • LLрасстояние до разлома, известное по геологии или сейсмике, м
  • k, ϕ, μ, ct, rwk,\ \phi,\ \mu,\ c_t,\ r_wпроницаемость (по удвоенной полке), пористость, вязкость, общая сжимаемость, радиус скважины, м², —, Па·с, 1/Па, м

Бірліктер: СИ; логарифм натуральный

Мына жағдайда орындалады
  • первая полка скрыта влиянием ствола, удвоенная полка сформировалась;
  • расстояние LL известно независимо — по такой записи его найти нельзя;
  • однородный пласт, разлом прямолинейный и непроницаемый.

Это вторая прямая предыдущего раздела, но без первой. Если о разломе не знать и обработать запись как обычный радиальный приток, наклон 2m2m примут за mm: проводимость получится вдвое меньше истинной, а скин — Sкаж=12[S+ln⁡(rw/L)]S_{\text{каж}} = \tfrac12\left[S + \ln(r_w/L)\right], то есть меньше истинного и часто отрицательным (вывод — в дополнительном разделе).

Басқа бірліктер жүйелерінде (1)
Метрические промысловые единицы
kh=1,866 qBμΔp2′,S=2,303[Δp1чm2−lg⁡kϕμctrw2+2,098+lg⁡2Lrw]kh = 1{,}866\,\frac{qB\mu}{\Delta p_2'},\qquad S = 2{,}303\left[\frac{\Delta p_{1\text{ч}}}{m_2} - \lg\frac{k}{\phi\mu c_tr_w^2} + 2{,}098 + \lg\frac{2L}{r_w}\right]

qq — м³/сут, μ\mu — мПа·с, kk — мД, ctc_t — 1/МПа, rwr_w и LL — м; m2m_2 — наклон прямой на десятичный цикл, Δp1ч\Delta p_{1\text{ч}} — снижение давления на ней через 1 ч, МПа.

1,866=2⋅0,9331{,}866 = 2\cdot0{,}933; множитель 2,3032{,}303 вместо обычного 1,1511{,}151 — потому что наклон удвоен, m2=2,303⋅Δp2′m_2 = 2{,}303\cdot\Delta p_2'; 2,0982{,}098 — как в формуле скина по кривой падения давления (КПД, Анализ кривой падения давления).

Диагностический график скважины в 10 метрах от разлома и производная той же скважины без разлома0,0010,010,11100,11Время после пуска t, чΔp и Δp′, МПаΔp2mₑmₑ: без разломаразлом в 10 м

График в двойном логарифмическом масштабе: время от 0,001 до 30 часов, давление от 0,02 до 3 МПа. Точки снижения давления круто растут до 0,5 МПа к 0,01 часа и дальше медленно до 1,3 МПа к 25 часам. Оранжевая производная проходит максимум около 0,2 МПа через 13 секунд, опускается до 0,078 МПа около 0,03 часа и затем поднимается к штриховой горизонтали 0,099 МПа, на которой лежит после 0,5 часа. Серая производная без разлома после того же горба спускается к нижней штриховой горизонтали 0,050 МПа.

ЕсептеуДиагностический график скважины в 10 м от разлома: снижение давления (точки) и производная по замерам (оранжевая линия); серая линия — производная той же скважины без разлома. После горба влияния ствола производная выходит сразу на удвоенный уровень 2me2m_e; первой полки mem_e в записи нет.
Есептеу мысалыСкважина у разлома, принятая за обычную
Бастапқы деректер (шартты)
Скважина
та же конструкция и тот же пласт, что в примерах выше, но в 10 м от разлома — по сейсмике
Запись
прямая от 2 до 20 ч: 1,0547 и 1,2829 МПа (ровно цикл); производная на полке около 0,0990 МПа
Шешім
  1. Обработка без учёта разлома: m=1,2829−1,0547=0,2282m = 1{,}2829 - 1{,}0547 = 0{,}2282 МПа на цикл; k=2,149⋅89,375/(0,2282⋅14)=60,1k = 2{,}149\cdot89{,}375/(0{,}2282\cdot14) = 60{,}1 мД.
  2. Δp1ч=1,0547−0,3⋅0,2282=0,9862\Delta p_{1\text{ч}} = 1{,}0547 - 0{,}3\cdot0{,}2282 = 0{,}9862 МПа; lg⁡kϕμctrw2=7,541\lg\dfrac{k}{\phi\mu c_tr_w^2} = 7{,}541; S=1,151 (0,9862/0,2282−7,541+2,098)=−1,29S = 1{,}151\,(0{,}9862/0{,}2282 - 7{,}541 + 2{,}098) = -1{,}29 — скважина выглядит обработанной.
  3. С учётом разлома наклон — это 2m2m: k=2⋅60,1=120,2k = 2\cdot60{,}1 = 120{,}2 мД. Скин: lg⁡kϕμctrw2=7,842\lg\dfrac{k}{\phi\mu c_tr_w^2} = 7{,}842, lg⁡(2L/rw)=lg⁡200=2,301\lg(2L/r_w) = \lg 200 = 2{,}301; S=2,303 (0,9862/0,2282−7,842+2,098+2,301)=2,02S = 2{,}303\,(0{,}9862/0{,}2282 - 7{,}842 + 2{,}098 + 2{,}301) = 2{,}02.
  4. Связь двух ответов: 12[S+ln⁡(rw/L)]=12 (2,02+ln⁡0,01)=−1,29\tfrac12[S + \ln(r_w/L)] = \tfrac12\,(2{,}02 + \ln 0{,}01) = -1{,}29 — ровно кажущийся скин из второго шага.
Жауап: Без знания о разломе получили бы k≈60,1k \approx 60{,}1 мД и скин −1,29; с учётом разлома — k≈120,2k \approx 120{,}2 мД и скин 2,02.

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

Восстановление давления у разлома

Чаще разломы ищут по кривой восстановления давления (КВД): дебит во время записи известен точно — он равен нулю. Расчёт требует теперь двух суперпозиций: во времени (пуск и остановка) и в пространстве (скважина и отражение). Для времени остановки Δt\Delta t после работы tpt_p:

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

Бірліктер: СИ; логарифм натуральный

Мына жағдайда орындалады
  • до пуска давление в пласте было pip_i; дебит постоянен в течение tpt_p;
  • для самой скважины — логарифмическое приближение, влияние ствола закончилось; других границ нет.

Первое слагаемое — обычная прямая Хорнера (Анализ кривой восстановления давления). Два слагаемых с E1E_1 — отражение, которое начало работать tp+Δtt_p + \Delta t назад и «остановилось» Δt\Delta t назад. Какое из них важнее, зависит от того, успело ли возмущение дойти до разлома за время работы, — от отношения L2/(ηtp)L^2/(\eta t_p).

Два случая КВД у разлома
Начало КВД
Разлом почувствовался при работе: L2<ηtpL^2 < \eta t_p
прямая Хорнера с верным наклоном mm, но p∗p^* ниже начального на meE1(L2/(ηtp))m_eE_1\big(L^2/(\eta t_p)\big)
Не почувствовался: L2>ηtpL^2 > \eta t_p
прямая Хорнера с p∗≈pip^* \approx p_i
Конец КВД
Разлом почувствовался при работе: L2<ηtpL^2 < \eta t_p
точки уходят вверх от прямой; при ηΔt≫L2\eta\Delta t \gg L^2 — наклон 2m2m
Не почувствовался: L2>ηtpL^2 > \eta t_p
небольшой провал под прямой и возврат к ней
Что даёт
Разлом почувствовался при работе: L2<ηtpL^2 < \eta t_p
расстояние по пересечению прямых Хорнера
Не почувствовался: L2>ηtpL^2 > \eta t_p
в лучшем случае признак границы и грубую оценку расстояния
Расстояние до разлома по пересечению прямых Хорнера
L=η Δtexγ,Δtex=tp Δtxtp+ΔtxL = \sqrt{\frac{\eta\,\Delta t_{ex}}{\gamma}},\qquad \Delta t_{ex} = \frac{t_p\,\Delta t_x}{t_p + \Delta t_x}
мұндағы:
  • Δtx\Delta t_xвремя остановки, при котором пересекаются продолжения первой и второй прямых Хорнера, с
  • Δtex\Delta t_{ex}соответствующее эквивалентное время (Влияние времени работы на анализ КВД), с
  • η, γ\eta,\ \gammaпьезопроводность и γ=1,781\gamma = 1{,}781, м²/с, —

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

Мына жағдайда орындалады
  • разлом почувствовался задолго до остановки: L2≪ηtpL^2 \ll \eta t_p, и E1(L2/(ηtp))E_1\big(L^2/(\eta t_p)\big) можно заменить логарифмом;
  • обе прямые Хорнера сформировались; недосформированная вторая прямая, как и при КПД, занижает LL.

Обозначим отношение Хорнера H=(tp+Δt)/ΔtH = (t_p + \Delta t)/\Delta t. Первая прямая — pi−pws=me(ln⁡H+ln⁡Hx)p_i - p_{ws} = m_e(\ln H + \ln H_x), где Hx=ηtp/(γL2)H_x = \eta t_p/(\gamma L^2); вторая — pi−pws=2meln⁡Hp_i - p_{ws} = 2m_e\ln H. Они пересекаются при H=HxH = H_x; выразив отсюда LL через Δtx\Delta t_x, получаем формулу. При Δtx≪tp\Delta t_x \ll t_p эквивалентное время почти равно Δtx\Delta t_x, и формула совпадает с формулой для КПД.

График Хорнера восстановления давления после пробной эксплуатации скважины у разлома100 00010 000100010010128,228,428,628,829,0Отношение Хорнера (tₚ + Δt)/ΔtЗабойное давление, МПаpᵢ = 29,0 МПапоздние точки выше прямой —сказывается разломp* = 28,88 МПа — ниже pᵢ

График Хорнера: по горизонтали отношение Хорнера от ста тысяч до единицы на логарифмической шкале, по вертикали давление от 28,2 до 29,1 МПа. Точки быстро поднимаются от 28,29 МПа, с отношения около 800 ложатся на штриховую прямую, а после отношения около 20 уходят вверх от неё и к отношению 2 доходят до 28,93 МПа. Прямая продолжена до отношения 1 и упирается в 28,88 МПа — ниже горизонтали начального давления 29 МПа.

ЕсептеуГрафик Хорнера КВД после 398 ч пробной эксплуатации (точки). Штриховая прямая проведена через замеры 1 и 5 ч и продолжена до отношения Хорнера 1, где читается p∗p^*; горизонталь вверху — начальное давление 29 МПа. Поздние точки уходят вверх от прямой.
Есептеу мысалыКВД после пробной эксплуатации: p∗p^* ниже начального давления
Бастапқы деректер (шартты)
Скважина
из примеров выше; после 398,1 ч работы её закрыли на устье
Замеры КВД
1 ч — 28,5831 МПа; 5,012 ч — 28,6626 МПа
Начальное давление
29,0 МПа (замерено до пуска)
Шешім
  1. Отношения Хорнера: (398,1+1)/1=399,1(398{,}1 + 1)/1 = 399{,}1 и (398,1+5,012)/5,012=80,43(398{,}1 + 5{,}012)/5{,}012 = 80{,}43; разность десятичных логарифмов 0,6957.
  2. Наклон: m=(28,6626−28,5831)/0,6957=0,1143m = (28{,}6626 - 28{,}5831)/0{,}6957 = 0{,}1143 МПа на цикл — почти тот же, что у первой прямой КПД: проницаемость по нему верная.
  3. Продолжение до отношения 1: p∗=28,6626+0,1143⋅lg⁡80,43=28,880p^* = 28{,}6626 + 0{,}1143\cdot\lg 80{,}43 = 28{,}880 МПа — на 0,120 МПа ниже начального.
  4. Ожидаемое занижение: me=0,1143/2,303=0,0496m_e = 0{,}1143/2{,}303 = 0{,}0496 МПа; при LL = 230 м и η\eta = 2447 м²/ч L2/(ηtp)=0,0543L^2/(\eta t_p) = 0{,}0543, E1(0,0543)=2,390E_1(0{,}0543) = 2{,}390, meE1=0,119m_eE_1 = 0{,}119 МПа.
Жауап: p∗p^* первой прямой — 28,880 МПа, на 0,120 МПа ниже начального давления 29,0 МПа; это вклад отражения, а не истощение пласта.

Бұл нені білдіреді. За 400 ч из пласта отобрано слишком мало, чтобы заметно снизить давление в залежи, а p∗p^* всё равно ниже начального: разлом удерживает воронку депрессии со своей стороны, а прямая Хорнера по ранним точкам об этом «не знает». Если бы начальное давление не было известно, занижение можно было бы принять за истощение или за признак связи с уже разрабатываемой частью залежи. Поздние точки этой КВД уходят вверх от прямой и к сотням часов приближаются к наклону 2m2m; по ним и по пересечению прямых Хорнера находят расстояние — с теми же оговорками о недосформированной второй прямой. Данные примера условные.

Если же за время работы возмущение до разлома не дошло, КВД видит его только по небольшому отклонению от прямой Хорнера. Отклонение растёт, проходит максимум примерно через Δt≈L2/η\Delta t \approx L^2/\eta после остановки и снова убывает. При ηtp≪L2\eta t_p \ll L^2 его наибольшее значение выражается просто:

Наибольшее отклонение КВД от прямой и предельное расстояние обнаружения разлома
δpmax⁡=qB tp4πe ϕcth L2⇒Lmax⁡=qB tp4πe ϕcth δp\delta p_{\max} = \frac{qB\,t_p}{4\pi\e\,\phi c_th\,L^2}\quad\Rightarrow\quad L_{\max} = \sqrt{\frac{qB\,t_p}{4\pi\e\,\phi c_th\,\delta p}}
мұндағы:
  • δpmax⁡\delta p_{\max}наибольшее отклонение давления КВД от прямой Хорнера бесконечного пласта, Па
  • δp\delta pнаименьшее отклонение, которое удаётся надёжно различить (разрешение и шум манометра), Па
  • q, B, tpq,\ B,\ t_pдебит перед остановкой, объёмный коэффициент, время работы, м³/с, м³/м³, с
  • ϕ, ct, h\phi,\ c_t,\ hпористость, общая сжимаемость, толщина пласта, —, 1/Па, м
  • e\eоснование натуральных логарифмов, 2,718, —

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

Мына жағдайда орындалады
  • разлом за время работы не почувствовался: ηtp≪L2\eta t_p \ll L^2 (практически годится уже при L2>3ηtpL^2 > 3\eta t_p);
  • запись длится не меньше L2/ηL^2/\eta после остановки;
  • однородный пласт; других помех — работы соседних скважин, дрейфа манометра — нет.

В формулу не входят ни проницаемость, ни вязкость: «сигнал» разлома определяется отобранным объёмом qBtpqBt_p, упругоёмкостью пласта ϕcth\phi c_th и расстоянием. От проницаемости зависит только, когда ждать максимума. Поэтому радиус, на котором КВД ещё замечает границу, задают не столько свойства пласта, сколько объём отбора перед остановкой и разрешение манометра.

Басқа бірліктер жүйелерінде (2)
Метрические промысловые единицы
δpmax⁡=qB tp820 ϕcth L2,Lmax⁡=0,0349qB tpϕcth δp\delta p_{\max} = \frac{qB\,t_p}{820\,\phi c_th\,L^2},\qquad L_{\max} = 0{,}0349\sqrt{\frac{qB\,t_p}{\phi c_th\,\delta p}}

qq — м³/сут, tpt_p — ч, ctc_t — 1/МПа, hh и LL — м; давления — МПа.

820=24⋅4πe820 = 24\cdot4\pi\e (24 — часов в сутках); 0,0349=1/8200{,}0349 = 1/\sqrt{820}.

Промысловые единицы (oilfield)
Lmax⁡=0,0827qB tpϕcth δpL_{\max} = 0{,}0827\sqrt{\frac{qB\,t_p}{\phi c_th\,\delta p}}

qq — STB/d, tpt_p — ч, ctc_t — 1/psi, hh и LL — ft, δp\delta p — psi.

0,0827=5,615/8200{,}0827 = \sqrt{5{,}615/820}, где 5,615 — кубических футов в барреле.

Есептеу мысалыУвидит ли короткое испытание разлом
Бастапқы деректер (шартты)
Испытание
приток 6 ч с дебитом 110 м³/сут, затем КВД; пласт тот же: BB = 1,25; ϕ\phi = 0,19; ct=1,4⋅10−3c_t = 1{,}4\cdot10^{-3} МПа⁻¹; hh = 14 м; η\eta ≈ 2447 м²/ч
Манометр
надёжно различает отклонения от 0,5 кПа (0,0005 МПа)
Разлом
в 230 м от скважины
Шешім
  1. ηtp=14 682\eta t_p = 14\,682 м², а L2=52 900L^2 = 52\,900 м²: за время притока разлом не почувствовался.
  2. Наибольшее отклонение: qBtp=110⋅1,25⋅6=825qBt_p = 110\cdot1{,}25\cdot6 = 825, ϕcth=0,19⋅1,4⋅10−3⋅14=3,724⋅10−3\phi c_th = 0{,}19\cdot1{,}4\cdot10^{-3}\cdot14 = 3{,}724\cdot10^{-3} м/МПа; δpmax⁡=825820⋅3,724⋅10−3⋅52 900=0,0051\delta p_{\max} = \dfrac{825}{820\cdot3{,}724\cdot10^{-3}\cdot52\,900} = 0{,}0051 МПа, то есть 5,1 кПа; наступает примерно через L2/η=22L^2/\eta = 22 ч после остановки.
  3. Расчёт по полной формуле КВД у разлома даёт 5,1 кПа через 19 ч — простая формула верна.
  4. Предельное расстояние при пороге 0,5 кПа: Lmax⁡=0,03498253,724⋅10−3⋅0,0005=735L_{\max} = 0{,}0349\sqrt{\dfrac{825}{3{,}724\cdot10^{-3}\cdot0{,}0005}} = 735 м; чтобы увидеть такой разлом, КВД должна длиться около Lmax⁡2/η=221L_{\max}^2/\eta = 221 ч.
Жауап: Разлом в 230 м даст отклонение около 5,1 кПа почти через сутки после остановки — его можно заметить, если КВД продлить; дальше 735 м это испытание разлом не покажет при любой длительности записи.

Бұл нені білдіреді. Отклонение в несколько килопаскалей — доли процента от роста давления при КВД; на графике Хорнера его видно только при хорошем манометре и спокойной записи. Чтобы «дотянуться» дальше, нужно больше отобрать: предельное расстояние растёт как корень из отобранного объёма, и вчетверо больший объём отбора удваивает радиус. Данные примера условные.

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

Пределы метода

  • Нужны обе прямые. Для метода пересечения запись должна дойти до ηt/L2≈10\eta t/L^2 \approx 10–25, то есть быть в десятки раз длиннее времени, когда разлом начал сказываться. Иначе вторая прямая положе настоящей и расстояние занижено (таблица выше). Ступень производной даёт расстояние раньше, но требует уверенности в модели.
  • Близкий разлом скрыт стволом. Если первая полка не успела проявиться, по записи нельзя отличить скважину у разлома от скважины в пласте вдвое меньшей проводимости; помогают только независимые сведения о разломе и сокращение влияния ствола.
  • Далёкий разлом не виден. При КПД — пока радиус исследования не приблизился к разлому; при КВД после короткой работы — дальше предельного расстояния обнаружения, сколько ни жди.
  • Удвоение — не доказательство. Ступень вдвое дают и другие модели: например, пласт, проводимость которого за некоторым радиусом вдвое меньше (переход похож, но не тождествен), или ложные «прямые» на участке влияния ствола. Проверка — геология, сейсмика, соседние скважины и согласие khkh и скина по обеим прямым.
  • Разлом может пропускать флюид. Тогда производная сначала поднимается, а затем возвращается к первой полке: через разлом начинается переток (Частично проводящие разломы).
  • Граница может быть не одна. Если производная поднимается выше удвоенного уровня или продолжает расти с наклоном ½ или 1, к первому разлому добавились другие границы (Несколько границ: параллельные разломы и клиновидные блоки).
  • Пласт за разломом и вдоль него. Метод предполагает однородный пласт и прямую границу через всю толщину. Изменение толщины или проницаемости вдоль границы, наклон плоскости разлома, различие проницаемости по направлениям меняют и время, и форму перехода; расстояние по простой модели — оценка, а не измерение.

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

Непроницаемый разлом
разлом, через который нет перетока флюида; для притока к скважине — прямолинейная непроницаемая граница.
Скважина-отражение
воображаемая скважина с той же историей дебита, поставленная зеркально за разломом на расстоянии 2L2L от реальной.
Удвоение наклона
переход от первой полулогарифмической прямой с наклоном mm ко второй с наклоном 2m2m; на производной — от полки mem_e к полке 2me2m_e.
Время пересечения txt_x
момент пересечения продолжений двух прямых; L=ktx/(γϕμct)L = \sqrt{kt_x/(\gamma\phi\mu c_t)}.
Ступень производной
подъём производной с mem_e до 2me2m_e по закону 1+e−L2/(ηt)1 + \e^{-L^2/(\eta t)}; занимает около 2,7 логарифмического цикла.
Полурадиальный приток
приток к скважине у самой непроницаемой границы; наклон прямой сразу удвоенный.
Предельное расстояние обнаружения
наибольшее расстояние до границы, при котором её влияние на КВД ещё превышает разрешение манометра; растёт как корень из отобранного объёма.

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

Қате

«Разлом проявится, как только до него дойдёт радиус исследования».

Шын мәнінде

На производной — да: когда радиус исследования равен LL, она уже на 2 % выше полки. В самом давлении вклад отражения заметен позже, когда радиус исследования приближается к 2L2L, а удвоение наклона устанавливается только при ηt/L2≈10\eta t/L^2 \approx 10 — в 40 раз позже.

Қате

«Вторая прямая — это последние точки записи».

Шын мәнінде

Только если к концу записи наклон действительно удвоился. В примере урока хорда двух последних замеров круче первой прямой лишь в 1,91 раза, и расстояние по ней занижено на 11 %.

Қате

«Отрицательный скин по исследованию означает, что скважина хорошо очищена».

Шын мәнінде

У скважины рядом с разломом, принятой за обычную, кажущийся скин равен 12[S−ln⁡(L/rw)]\tfrac12[S - \ln(L/r_w)], а проницаемость занижена вдвое. Отрицательный скин у необработанной скважины — повод проверить геологию.

Қате

«Чем длиннее КВД, тем дальше она «видит»».

Шын мәнінде

Только до предела: наибольшее отклонение от прямой Хорнера определяется объёмом, отобранным перед остановкой, и расстоянием. Если оно меньше разрешения манометра, никакая длительность КВД разлом не покажет.

Қате

«p∗p^* по прямой Хорнера в новой залежи равно начальному давлению».

Шын мәнінде

Если разлом почувствовался за время работы, p∗p^* первой прямой ниже начального на meE1(L2/(ηtp))m_eE_1\big(L^2/(\eta t_p)\big) — в примере урока на 0,12 МПа.

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

Выводы: ступень производной, кажущийся скин, отклонение КВД

Ступень производной. Производная интегральной показательной функции  dE1(x)/ dx=−e−x/x\dd E_1(x)/\dd x = -\e^{-x}/x. Для x=a/tx = a/t имеем  dx/ dln⁡t=t  dx/ dt=−x\dd x/\dd\ln t = t\,\dd x/\dd t = -x, поэтому  dE1/ dln⁡t\dd E_1/\dd\ln t равна (−e−x/x)⋅(−x)=e−x(-\e^{-x}/x)\cdot(-x) = \e^{-x}, то есть e−a/t\e^{-a/t}. Вклад отражения в давление — meE1(L2/(ηt))m_eE_1\big(L^2/(\eta t)\big), поэтому его вклад в производную — mee−L2/(ηt)m_e\e^{-L^2/(\eta t)}, а скважина после окончания влияния ствола даёт mem_e.

Кажущийся скин. Обозначим ℓ=ln⁡(2L/rw)\ell = \ln(2L/r_w). На удвоенной прямой Δp=2me(ln⁡tD+0,809+S−ℓ)\Delta p = 2m_e(\ln t_D + 0{,}809 + S - \ell). Если её принять за обычный радиальный приток Δp=me′(ln⁡tD′+0,809+2S′)\Delta p = m_e'(\ln t_D' + 0{,}809 + 2S'), то me′=2mem_e' = 2m_e и проницаемость вдвое меньше, а значит, tD′=tD/2t_D' = t_D/2. Приравнивая, получаем 2S′=S−ℓ+ln⁡22S' = S - \ell + \ln 2, а это S−ln⁡(L/rw)S - \ln(L/r_w), то есть S′=12[S+ln⁡(rw/L)]S' = \tfrac12\big[S + \ln(r_w/L)\big]. Для скважины радиусом 0,1 м в 10 м от разлома ln⁡(rw/L)=−4,6\ln(r_w/L) = -4{,}6: даже скин +2 превращается в −1,3.

Наибольшее отклонение КВД. Обозначим a=L2/ηa = L^2/\eta и пусть ηtp≪L2\eta t_p \ll L^2. При Δt≫tp\Delta t \gg t_p разность E1(a/(tp+Δt))−E1(a/Δt)E_1\big(a/(t_p + \Delta t)\big) - E_1(a/\Delta t) — это приращение E1(a/t)E_1(a/t) за время tpt_p, то есть примерно tp e−a/Δt/Δtt_p\,\e^{-a/\Delta t}/\Delta t. Эта величина наибольшая при Δt=a\Delta t = a и равна tp/(e a)t_p/(\e\,a). Умножив на me=qBμ/(4πkh)m_e = qB\mu/(4\pi kh) и подставив a=ϕμctL2/ka = \phi\mu c_tL^2/k, получаем δpmax⁡=qBtp/(4πe ϕcthL2)\delta p_{\max} = qBt_p/(4\pi\e\,\phi c_thL^2): вязкость и проницаемость сократились.

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

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

  1. 1

    Как меняются наклон полулогарифмической прямой и уровень логарифмической производной, когда возмущение давления полностью «почувствовало» один прямолинейный непроницаемый разлом?

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

    Дұрыс жауап: Б. Оба удваиваются

    Вклад отражения при ηt≫L2\eta t \gg L^2 сам растёт как meln⁡tm_e\ln t, то есть добавляет к наклону ещё mem_e на единицу натурального логарифма. Производная — это тот же наклон, поэтому и она выходит на 2me2m_e. Падение до нуля — признак границы постоянного давления.

  2. 2

    По первой прямой kk = 60 мД; ϕ\phi = 0,21, μ\mu = 1,1 мПа·с, ct=1,6⋅10−3c_t = 1{,}6\cdot10^{-3} МПа⁻¹. Продолжения первой и второй прямых пересеклись через 100 ч после пуска. Найдите расстояние до разлома.

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

    ϕμct=0,21⋅1,1⋅1,6⋅10−3=3,70⋅10−4\phi\mu c_t = 0{,}21\cdot1{,}1\cdot1{,}6\cdot10^{-3} = 3{,}70\cdot10^{-4}. L=0,044760⋅100/(3,70⋅10−4)=0,04471,62⋅107≈0,0447⋅4030≈180L = 0{,}0447\sqrt{60\cdot100/(3{,}70\cdot10^{-4})} = 0{,}0447\sqrt{1{,}62\cdot10^7} \approx 0{,}0447\cdot4030 \approx 180 м. Если вторая прямая проведена по ещё не удвоившемуся наклону, настоящее расстояние больше.

  3. 3

    Уровень первой полки производной — 0,040 МПа, через 30 ч после пуска производная равна 0,062 МПа, η\eta = 1500 м²/ч. Найдите расстояние до разлома и оцените, когда производная дойдёт до 0,076 МПа (90 % ступени).

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

    Δp′/me=1,55\Delta p'/m_e = 1{,}55; ln⁡[1/(1,55−1)]=ln⁡1,818=0,598\ln[1/(1{,}55 - 1)] = \ln 1{,}818 = 0{,}598; L=1500⋅30⋅0,598≈26 900≈164L = \sqrt{1500\cdot30\cdot0{,}598} \approx \sqrt{26\,900} \approx 164 м. Отношение 1,9 наступает при ηt/L2=1/ln⁡(1/0,9)=9,5\eta t/L^2 = 1/\ln(1/0{,}9) = 9{,}5, то есть через 9,5⋅26 900/1500≈1709{,}5\cdot26\,900/1500 \approx 170 ч.

  4. 4

    Скважина радиусом 0,1 м стоит в 8 м от разлома, влияние ствола длится около 2 ч. Обработка КПД как обычного радиального притока дала kk = 40 мД и скин −1,8. Что вероятнее всего на самом деле?

    1. А.Запись обработать нельзя
    2. Б.Проницаемость около 80 мД, скин около +0,8
    3. В.Проницаемость около 20 мД, скин около −3,6
    4. Г.Скважина действительно стимулирована, kk = 40 мД
    Жауабын көрсету

    Дұрыс жауап: Б. Проницаемость около 80 мД, скин около +0,8

    Первая полка скрыта стволом, видна только удвоенная: наклон 2m2m приняли за mm, поэтому истинная проницаемость вдвое больше. Кажущийся скин S′=12[S+ln⁡(rw/L)]S' = \tfrac12[S + \ln(r_w/L)], откуда S=2S′−ln⁡(rw/L)=−3,6+ln⁡80≈−3,6+4,38≈0,8S = 2S' - \ln(r_w/L) = -3{,}6 + \ln 80 \approx -3{,}6 + 4{,}38 \approx 0{,}8.

  5. 5

    Во сколько раз нужно увеличить объём, отобранный перед КВД, чтобы предельное расстояние обнаружения разлома выросло вдвое? Изменится ли время, через которое ждать наибольшего отклонения от прямой Хорнера для разлома на этом новом расстоянии?

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

    Lmax⁡L_{\max} пропорционально корню из qBtpqBt_p, поэтому объём нужен вчетверо больший. Максимум отклонения наступает примерно через L2/ηL^2/\eta после остановки: для вдвое более далёкого разлома — вчетверо позже. Для того же разлома время максимума от объёма отбора не зависит, растёт только величина отклонения.

  6. 6

    Производная поднялась над первой полкой, но через полтора цикла вернулась к прежнему уровню. Что это скорее всего?

    1. А.Граница постоянного давления
    2. Б.Замкнутая область
    3. В.Частично проводящий разлом
    4. Г.Непроницаемый разлом
    Жауабын көрсету

    Дұрыс жауап: В. Частично проводящий разлом

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

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