Skip to content
Brent$117,6Газ TTF$1 192USD453,87 ₸
MunaiHub
Физика пластаРежимы залежей и материальный баланс Продвинутый 21 мин чтения

Линейная форма материального баланса (метод Хавлены–Одеха)

Метод Хавлены–Одеха (Havlena–Odeh method), или линейная форма материального баланса, — запись уравнения материального баланса нефтяной залежи как уравнения прямой. Каждая дата замера пластового давления даёт на графике точку, координаты которой вычисляют по добыче и PVT-свойствам флюидов. Если представление о залежи верно, точки ложатся на прямую; её наклон и отрезок на оси дают геологические запасы нефти, размер газовой шапки или параметр водоносного горизонта, а отклонение точек от прямой показывает, что в представлении о залежи чего-то не хватает.

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

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

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

  • Сокращённая запись баланса: F=N(Eo+mEg+Efw)+WeF = N(E_o + mE_g + E_{fw}) + W_e. Слева — отбор, вычисленный по добыче, справа — неизвестные NN, mm и WeW_e, умноженные на расширения, вычисленные по PVT-свойствам при замеренном давлении.
  • Каждая дата замера давления — точка на графике. Координаты выбирают так, чтобы при верной модели залежи точки легли на прямую: её наклон и отрезок — искомые величины, а систематическое отклонение точек от прямой — признак, что модель неполна или данные ошибочны.
  • Залежь без газовой шапки: кажущиеся запасы F/(Eo+Efw)F/(E_o + E_{fw}) постоянны, если притока воды нет, и растут со временем, если вода входит. Постоянство ещё не доказывает замкнутость: неучтённая шапка или небольшой водоносный горизонт выше давления насыщения дают постоянные, но завышенные запасы.
  • С газовой шапкой неизвестного размера график F/EoF/E_o от Eg/EoE_g/E_o даёт отрезок NN и наклон mNmN. Но EoE_o и EgE_g растут с падением давления почти пропорционально, точки собираются в узкую полосу, и прямая плохо разделяет NN и mm: размер шапки обычно берут по геологическим данным, а график служит проверкой.
  • С водоносным горизонтом приток выражают через модель горизонта: WeW_e равен её параметру, умноженному на функцию истории давления. Тогда F/EtF/E_t от этой функции, делённой на EtE_t, — прямая с отрезком NN и наклоном, равным параметру горизонта, если модель выбрана верно; изгиб графика означает неверную модель, а не неточный коэффициент.
  • Индексы режимов — доли пластового объёма отобранных нефти и газа, которые дали расширение нефти, газовая шапка, вода и поры и приток воды. Их сумма равна единице, когда точка лежит на прямой; отличие суммы от единицы — относительное отклонение точки.
Ключевые соотношения
F=N (Eo+m Eg+Efw)+WeF = N\,(E_o + m\,E_g + E_{fw}) + W_e

сокращённая запись уравнения материального баланса нефтяной залежи

FEo+Efw=N+WeEo+Efw\frac{F}{E_o + E_{fw}} = N + \frac{W_e}{E_o + E_{fw}}

кажущиеся запасы залежи без газовой шапки: постоянны, если притока воды нет

FEo+Efw0=N+mN Eg+Efw0Eo+Efw0\frac{F}{E_o + E_{fw0}} = N + mN\,\frac{E_g + E_{fw0}}{E_o + E_{fw0}}

залежь с газовой шапкой без притока воды; Efw=(1+m) Efw0E_{fw} = (1 + m)\,E_{fw0}

FEt=N+Ja 1Et∫0t(pi−p) dt\frac{F}{E_t} = N + J_a\,\frac{1}{E_t}\int_0^t (p_i - p)\dd t

залежь со стационарным водоносным горизонтом; Et=Eo+mEg+EfwE_t = E_o + mE_g + E_{fw}

NEoH+NmEgH+NEfwH+We−WpBwH=1\frac{NE_o}{H} + \frac{NmE_g}{H} + \frac{NE_{fw}}{H} + \frac{W_e - W_pB_w}{H} = 1

индексы режимов; HH — пластовый объём отобранных нефти и газа

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

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

  • записать уравнение материального баланса в сокращённых обозначениях FF, EoE_o, EgE_g, EfwE_{fw}, WeW_e и привести его к уравнению прямой для залежи без газовой шапки, с газовой шапкой и с водоносным горизонтом;
  • по графику кажущихся запасов от накопленной добычи отличить замкнутую залежь от залежи с притоком воды и распознать неучтённую газовую шапку;
  • определить запасы нефти и размер газовой шапки по графику и объяснить, почему их совместное определение требует очень точных замеров давления;
  • определить запасы и коэффициент продуктивности водоносного горизонта по прямой для стационарной модели притока и распознать неверную модель горизонта по форме графика;
  • рассчитать индексы режимов по истории разработки и объяснить, почему их сумма равна единице.

Сокращённая запись и идея прямой

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

Уравнение материального баланса в сокращённой записи
F=N (Eo+m Eg+Efw)+WeF = N\,(E_o + m\,E_g + E_{fw}) + W_e
где:
  • FFпластовый объём отбора нефти, газа и воды при текущем среднем давлении (underground withdrawal), м³
  • NNгеологические запасы нефти в стандартных условиях, м³
  • EoE_oрасширение нефти с растворённым в ней газом на 1 м³ запасов, м³/м³
  • EgE_gрасширение газовой шапки на 1 м³ запасов нефти при m=1m = 1, м³/м³
  • EfwE_{fw}расширение связанной воды и уменьшение объёма пор на 1 м³ запасов, м³/м³
  • mmотносительный размер газовой шапки, —
  • WeW_eнакопленный приток законтурной воды в пластовых условиях, м³

Единицы: согласованные единицы: объёмы в м³ (запасы и добыча — в стандартных условиях, 0,101325 МПа и 20 °C, остальное — в пластовых); в промысловых единицах — STB, scf и rb

Справедливо, если
  • допущения общего уравнения: залежь — один «бак» с единым средним давлением, температура постоянна, свойства флюидов описывает модель «чёрной нефти»;
  • закачки нет; закачанные объёмы WinjBw+GinjBgW_{inj}B_w + G_{inj}B_g вычитают из FF.

Слева — то, что вычисляют по замерам добычи и давления. Справа — неизвестные (NN, mm, WeW_e), умноженные на величины, которые вычисляют по PVT-свойствам при том же давлении. Метод состоит в том, чтобы разделить уравнение на подходящую группу и получить прямую, у которой неизвестные стоят в наклоне и в отрезке.

Слагаемые сокращённой записи
F=Np[Bo+(Rp−Rs)Bg]+WpBwEo=(Bo−Boi)+(Rsi−Rs) BgEg=Boi(BgBgi−1)Efw=(1+m) Efw0,Efw0=Boi cwSw+cf1−Sw Δp\begin{aligned} F &= N_p\big[B_o + (R_p - R_s)B_g\big] + W_pB_w\\ E_o &= (B_o - B_{oi}) + (R_{si} - R_s)\,B_g\\ E_g &= B_{oi}\left(\frac{B_g}{B_{gi}} - 1\right)\\ E_{fw} &= (1 + m)\,E_{fw0},\qquad E_{fw0} = B_{oi}\,\frac{c_wS_w + c_f}{1 - S_w}\,\Delta p\end{aligned}
где:
  • Np, WpN_p,\ W_pнакопленная добыча нефти и воды в стандартных условиях, м³
  • RpR_pнакопленный газовый фактор: весь добытый газ на всю добытую нефть, м³/м³
  • Bo, Rs, Bg, BwB_o,\ R_s,\ B_g,\ B_wобъёмные коэффициенты нефти, газа и воды и газосодержание при текущем среднем давлении, м³/м³
  • Boi, Rsi, BgiB_{oi},\ R_{si},\ B_{gi}те же свойства при начальном давлении, м³/м³
  • Efw0E_{fw0}расширение связанной воды и уменьшение объёма пор при m=0m = 0, то есть для одной нефтяной части, м³/м³
  • cw, cf, Swc_w,\ c_f,\ S_wсжимаемость воды и порового объёма, насыщенность связанной водой, 1/МПа, 1/МПа, —
  • Δp\Delta pснижение среднего пластового давления от начального, pi−pp_i - p, МПа

Единицы: FF — м³; расширения — м³ на 1 м³ запасов нефти; давление и сжимаемости — МПа и 1/МПа (или psi и 1/psi)

Справедливо, если
  • все свойства взяты при одном и том же среднем давлении pp;
  • насыщенность связанной водой и сжимаемость пор в газовой шапке такие же, как в нефтяной части.

Вывод каждого слагаемого — в уроке Общее уравнение материального баланса нефтяной залежи. Новое здесь одно: слагаемое воды и пор пропорционально объёму пор всей залежи, поэтому оно делится между нефтяной частью (Efw0E_{fw0}) и шапкой (mEfw0mE_{fw0}). Отсюда точная группировка F=N(Eo+Efw0)+mN(Eg+Efw0)+WeF = N(E_o + E_{fw0}) + mN(E_g + E_{fw0}) + W_e: каждая зона — со своими водой и порами.

Уравнение связывает четыре группы величин: отбор, расширения, неизвестные запасы с размером шапки и приток воды. Для одной даты это одно уравнение, и найти из него можно только одно неизвестное — при условии, что остальные известны. Но замеров давления за годы разработки много, и каждый даёт своё уравнение с теми же NN и mm. Если поделить уравнение на подходящую группу, оно превращается в y=a+bxy = a + bx, где yy и xx вычисляются для каждой даты, а неизвестные стоят в отрезке aa и наклоне bb.

Какие координаты выбирают
Залежь
без шапки и без воды
Ось yy
FF
Ось xx
Eo+EfwE_o + E_{fw}
Наклон
NN
Отрезок
0
Залежь
без шапки, проверка на приток
Ось yy
F/(Eo+Efw)F/(E_o + E_{fw})
Ось xx
NpN_p или время
Наклон
—
Отрезок
NN
Залежь
шапка известного размера
Ось yy
FF
Ось xx
Eo+mEg+EfwE_o + mE_g + E_{fw}
Наклон
NN
Отрезок
0
Залежь
шапка неизвестного размера
Ось yy
F/EoF/E_o
Ось xx
Eg/EoE_g/E_o
Наклон
mNmN
Отрезок
NN
Залежь
водоносный горизонт
Ось yy
F/EtF/E_t
Ось xx
функция притока, делённая на EtE_t
Наклон
параметр горизонта
Отрезок
NN

Et=Eo+mEg+EfwE_t = E_o + mE_g + E_{fw} — полное расширение на 1 м³ запасов. В строке с шапкой неизвестного размера связанная вода и поры опущены; как их учесть, показано в разделе о газовой шапке.

Порядок работы
  1. 1
    Собрать данные на даты замеров

    Среднее пластовое давление, приведённое к одной глубине, накопленные добычу нефти, газа и воды и закачку.

  2. 2
    Найти свойства флюидов

    BoB_o, RsR_s, BgB_g при каждом давлении — по PVT-исследованиям, приведённым к условиям сепарации (Модель «чёрной нефти»: газосодержание и объёмные коэффициенты).

  3. 3
    Вычислить отбор и расширения

    FF, EoE_o, EgE_g и EfwE_{fw} для каждой даты.

  4. 4
    Выбрать координаты

    По предполагаемым источникам энергии — по таблице выше.

  5. 5
    Построить график и оценить форму

    Прямая — модель согласуется с историей, наклон и отрезок дают параметры. Изгиб или разброс — модель неполна или данные ошибочны.

  6. 6
    Проверить вывод

    Сравнить запасы с объёмным подсчётом, размер шапки — с геологическими данными, приток воды — с обводнением скважин; рассчитать индексы режимов.

Залежь без газовой шапки: есть ли приток воды

Если газовой шапки нет (m=0m = 0, и тогда Efw=Efw0E_{fw} = E_{fw0}), уравнение принимает вид F=N(Eo+Efw)+WeF = N(E_o + E_{fw}) + W_e. Для замкнутой залежи We=0W_e = 0, и точки «FF — Eo+EfwE_o + E_{fw}» ложатся на прямую, проходящую через начало координат, с наклоном NN (пример с двумя датами — в уроке о материальном балансе). Удобнее, однако, разделить обе части на расширение.

Кажущиеся запасы нефти залежи без газовой шапки
FEo+Efw=N+WeEo+Efw\frac{F}{E_o + E_{fw}} = N + \frac{W_e}{E_o + E_{fw}}
где:
  • F/(Eo+Efw)F/(E_o + E_{fw})кажущиеся запасы: запасы, вычисленные по одной дате в предположении, что притока воды нет, м³
  • NNгеологические запасы нефти, м³
  • WeW_eнакопленный приток воды к той же дате, м³
  • Eo, EfwE_o,\ E_{fw}расширение нефти с её газом и связанной воды с порами на 1 м³ запасов, м³/м³

Единицы: м³ (в промысловых единицах — STB)

Справедливо, если
  • газовой шапки нет; закачки нет или она вычтена из FF;
  • давление снизилось настолько, что его погрешность мала по сравнению с Δp\Delta p.

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

  • Горизонталь на уровне объёмной оценки — залежь замкнута, и вся её нефть сообщается со скважинами.
  • Подъём — приток законтурной воды. У большого водоносного горизонта кривая продолжает расти. У небольшого рост замедляется: упругий запас горизонта исчерпывается, а ниже давления насыщения расширение нефти с газом растёт быстрее притока, и кривая может даже пойти вниз (Модели притока законтурной воды).
  • Почти горизонталь заметно выше объёмной оценки — есть источник энергии, который меняется с давлением примерно так же, как расширение самой нефти: неучтённая газовая шапка (у нефти урока такая кривая медленно снижается), а выше давления насыщения — и небольшой водоносный горизонт. Либо занижен объёмный подсчёт.
  • Устойчивое снижение около объёмной оценки и ниже источником энергии не объясняется: обычно это ошибки замеров давления или учёта добычи либо отток флюидов в соседние объекты.
Кажущиеся запасы в долях объёмной оценки в зависимости от доли отобранных запасов: замкнутая залежь, залежь с водоносным горизонтом и залежь с неучтённой газовой шапкой05101520024681012Отобрано, % объёмных запасовКажущиеся запасы, доли объёмныхзамкнутая залежьгазовая шапка не учтенаводоносный горизонт

График: по горизонтали доля отобранных запасов от 0 до 20 процентов, по вертикали кажущиеся запасы в долях объёмной оценки от 0 до 12. Горизонтальная линия на уровне 1 с двумя точками — замкнутая залежь. Кривая залежи с водоносным горизонтом поднимается от 2,9 при 2 процентах до 10,5 при 20 процентах. Кривая залежи с неучтённой газовой шапкой идёт почти горизонтально, медленно снижаясь от 2,9 до 2,6.

РасчётКажущиеся запасы F/(Eo+Efw)F/(E_o + E_{fw}) в долях объёмной оценки в зависимости от доли отобранных запасов для трёх условных залежей урока: замкнутой (точки — примеры урока о материальном балансе), с водоносным горизонтом и с газовой шапкой, которую при расчёте «забыли». Приток воды поднимает кривую, неучтённая шапка держит её почти горизонтально, но в 2,6–2,9 раза выше единицы.

Залежь с газовой шапкой

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

Залежь с газовой шапкой известного размера
F=N [(Eo+Efw0)+m (Eg+Efw0)]F = N\,\big[(E_o + E_{fw0}) + m\,(E_g + E_{fw0})\big]
где:
  • FFпластовый объём отбора, м³
  • NNгеологические запасы нефти, м³
  • mmотносительный размер шапки — известен по геологическим данным, —
  • Eo, Eg, Efw0E_o,\ E_g,\ E_{fw0}расширение нефти с её газом, газа шапки (при m=1m = 1) и связанной воды с порами (при m=0m = 0) на 1 м³ запасов, м³/м³

Единицы: м³

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

При верном mm точки «FF — квадратная скобка» ложатся на прямую через начало координат с наклоном NN, а запасы, вычисленные по каждой дате, от даты не зависят. При неверном mm точки систематически отклоняются от прямой: вычисленные по датам запасы дрейфуют. Если пренебречь водой и порами, получится запись F=N(Eo+mEg)F = N(E_o + mE_g).

Размер шапки определяют по геологическим данным: по положению газонефтяного и водонефтяного контактов, картам эффективных толщин и каротажу. mm — отношение начального пластового объёма газа шапки к начальному пластовому объёму нефти; в стандартных кубометрах в шапке G=mNBoi/BgiG = mNB_{oi}/B_{gi} газа. Если надёжных геологических данных нет, оба неизвестных можно попытаться найти из самого баланса, разделив обе его части на Eo+Efw0E_o + E_{fw0}.

Залежь с газовой шапкой неизвестного размера
FEo+Efw0=N+mN Eg+Efw0Eo+Efw0\frac{F}{E_o + E_{fw0}} = N + mN\,\frac{E_g + E_{fw0}}{E_o + E_{fw0}}
где:
  • F/(Eo+Efw0)F/(E_o + E_{fw0})ордината точки, м³
  • (Eg+Efw0)/(Eo+Efw0)(E_g + E_{fw0})/(E_o + E_{fw0})абсцисса точки, —
  • NNотрезок на оси ординат — геологические запасы нефти, м³
  • mNmNнаклон — начальный пластовый объём шапки, делённый на BoiB_{oi}, м³

Единицы: м³ (в промысловых единицах — STB)

Справедливо, если
  • притока воды нет;
  • снижение давления велико по сравнению с погрешностью его замера.

Если пренебречь водой и порами, это знакомая запись «F/EoF/E_o от Eg/EoE_g/E_o». Отрезок — запасы нефти, наклон, делённый на отрезок, — размер шапки mm. Если точки не ложатся на прямую, в залежь, скорее всего, входит вода.

В других системах единиц (1)
Через газ шапки в стандартных условиях
FEo=N+G Bg−BgiEo\frac{F}{E_o} = N + G\,\frac{B_g - B_{gi}}{E_o}

G=mNBoi/BgiG = mNB_{oi}/B_{gi} — газ шапки, м³ в стандартных условиях; Bg−BgiB_g - B_{gi} — расширение 1 м³ этого газа, м³/м³; связанная вода и поры опущены.

Постоянных нет: mNEg=mNBoi(Bg/Bgi−1)=G (Bg−Bgi)mNE_g = mNB_{oi}(B_g/B_{gi} - 1) = G\,(B_g - B_{gi}). Наклон в этой записи — GG, а не mNmN: важно не смешивать два определения расширения газа.

История разработки условной залежи с газовой шапкой
Год
1-й год
pp, МПа
17,47
NpN_p, тыс. м³
150
RpR_p, м³/м³
112
Год
2-й год
pp, МПа
16,81
NpN_p, тыс. м³
340
RpR_p, м³/м³
118
Год
3-й год
pp, МПа
16,04
NpN_p, тыс. м³
560
RpR_p, м³/м³
128
Год
4-й год
pp, МПа
15,21
NpN_p, тыс. м³
790
RpR_p, м³/м³
140
Год
5-й год
pp, МПа
14,32
NpN_p, тыс. м³
1020
RpR_p, м³/м³
155
Год
6-й год
pp, МПа
13,44
NpN_p, тыс. м³
1240
RpR_p, м³/м³
170

Та же нефть, что в уроке Модель «чёрной нефти»: газосодержание и объёмные коэффициенты: давление насыщения 18 МПа, начальное газосодержание 110 м³/м³. Начальное давление на газонефтяном контакте равно давлению насыщения; притока воды нет, воду не добывают. RpR_p — накопленный газовый фактор.

Расчётный примерЗапасы нефти и размер газовой шапки
Исходные данные (условные)
История
таблица выше: давление, накопленная добыча нефти и накопленный газовый фактор по годам
Начальные свойства
Boi=1,35B_{oi} = 1{,}35 м³/м³; Rsi=110R_{si} = 110 м³/м³; Bgi=0,005983B_{gi} = 0{,}005983 м³/м³
Вода и поры
Sw=0,25S_w = 0{,}25; cw=0,45⋅10−3c_w = 0{,}45\cdot10^{-3} и cf=0,8⋅10−3c_f = 0{,}8\cdot10^{-3} МПа⁻¹, то есть (cwSw+cf)/(1−Sw)=1,2167⋅10−3(c_wS_w + c_f)/(1 - S_w) = 1{,}2167\cdot10^{-3} МПа⁻¹
Геологические данные
объёмный подсчёт нефтяной части — около 6,5 млн м³; размер шапки по контактам и картам m≈0,5m \approx 0{,}5
Решение
  1. Шестой год, p=13,44p = 13{,}44 МПа. Свойства при этом давлении: Rs=84,57R_s = 84{,}57 м³/м³, Bo=1,2806B_o = 1{,}2806 м³/м³, Bg=0,008013B_g = 0{,}008013 м³/м³. Отбор: F=1240⋅[1,2806+(170−84,57)⋅0,008013]=2437F = 1240\cdot\big[1{,}2806 + (170 - 84{,}57)\cdot0{,}008013\big] = 2437 тыс. м³.
  2. Расширения на кубометр запасов: Eo=(1,2806−1,35)+(110−84,57)⋅0,008013=0,1344E_o = (1{,}2806 - 1{,}35) + (110 - 84{,}57)\cdot0{,}008013 = 0{,}1344; Eg=1,35⋅(0,008013/0,005983−1)=0,4580E_g = 1{,}35\cdot(0{,}008013/0{,}005983 - 1) = 0{,}4580; Efw0=1,35⋅1,2167⋅10−3⋅(18−13,44)=0,00749E_{fw0} = 1{,}35\cdot1{,}2167\cdot10^{-3}\cdot(18 - 13{,}44) = 0{,}00749 м³/м³.
  3. Если считать, что шапки нет: F/(Eo+Efw0)=17,2F/(E_o + E_{fw0}) = 17{,}2 млн м³ — в 2,6 раза больше объёмной оценки нефтяной части. Без газовой шапки историю не объяснить.
  4. С шапкой по геологическим данным, m=0,5m = 0{,}5: N=F/[(Eo+Efw0)+0,5 (Eg+Efw0)]=6,50N = F/\big[(E_o + E_{fw0}) + 0{,}5\,(E_g + E_{fw0})\big] = 6{,}50 млн м³. По шести датам — от 6,47 до 6,50 млн м³ (таблица после примера): от даты запасы почти не зависят и согласуются с объёмной оценкой.
  5. Совместный подбор без геологических данных. Для каждой даты вычисляем x=(Eg+Efw0)/(Eo+Efw0)x = (E_g + E_{fw0})/(E_o + E_{fw0}) и y=F/(Eo+Efw0)y = F/(E_o + E_{fw0}); для шестого года x=3,280x = 3{,}280, y=17,17y = 17{,}17 млн м³. Прямая по шести точкам (метод наименьших квадратов): отрезок N=7,1N = 7{,}1 млн м³, наклон mN=3,07mN = 3{,}07 млн м³, то есть m=0,43m = 0{,}43.
  6. Качество двух описаний: точки отклоняются от подобранной прямой не больше чем на 0,16 %, от прямой с N=6,5N = 6{,}5 млн м³ и m=0,5m = 0{,}5 — не больше чем на 0,5 %. Оба описания укладываются в точность исходных данных.
  7. Чувствительность: если бы давления были известны лишь с точностью до 0,1 МПа (17,5; 16,8; 16,0; 15,2; 14,3; 13,4 МПа), тот же подбор дал бы N=1,6N = 1{,}6 млн м³ и m=2,9m = 2{,}9. Сильнее всего результат сдвигает первая точка: снижение давления 0,53 МПа превратилось в 0,5 МПа — ошибка 6 %, а её координаты делятся на самое малое расширение. Но и без неё подбор даёт N=7,6N = 7{,}6 млн м³ и m=0,38m = 0{,}38.
Ответ: При размере шапки 0,5 по геологическим данным запасы нефти — около 6,5 млн м³, и история это подтверждает. Без геологических данных баланс одинаково хорошо описывают пара 6,5 млн м³ и 0,5 и пара 7,1 млн м³ и 0,43.

Что это значит. Абсциссы шести точек меняются всего от 3,835 до 3,280: расширения нефти и газа шапки растут почти пропорционально. Отрезок на оси ординат находят, продолжая прямую примерно в 6 раз дальше ширины облака точек, и малое изменение наклона заметно сдвигает его. Надёжно определена средняя ордината облака — N(1+mxˉ)≈18,0N(1 + m\bar x) \approx 18{,}0 млн м³: обе пары дают её с точностью до долей процента. Поэтому размер шапки берут по геологическим данным, а баланс проверяет их согласие с историей; разделить NN и mm по одному балансу можно только при длинной истории, большом снижении давления и очень точных его замерах. Данные примера условные: история рассчитана по уравнению баланса для 6,5 млн м³ нефти и m=0,5m = 0{,}5, давления округлены до 0,01 МПа.

Замеры залежи с газовой шапкой в координатах графика для шапки неизвестного размера и две прямые: подобранная по точкам и отвечающая геологическим данным0123448121620(Eg + Efw0) / (Eo + Efw0)F / (Eo + Efw0), млн м³подбор по точкам: N = 7,1по геологии, m = 0,5: N = 6,5шесть замеров

График: по горизонтали абсцисса от 0 до 4, по вертикали ордината от 4 до 20 млн кубометров. Шесть точек собраны у правого края при абсциссах от 3,3 до 3,8 и ординатах от 17,2 до 18,9. Через них проходят две почти совпадающие прямые; продолженные штрихом влево, они расходятся и пересекают ось ординат в точках 6,5 и 7,1 млн кубометров.

РасчётШесть замеров залежи с газовой шапкой в координатах x=(Eg+Efw0)/(Eo+Efw0)x = (E_g + E_{fw0})/(E_o + E_{fw0}), y=F/(Eo+Efw0)y = F/(E_o + E_{fw0}) и две прямые, которые описывают их с точностью до долей процента: подобранная по точкам и отвечающая геологическим данным (m=0,5m = 0{,}5). Точки занимают узкую полосу абсцисс, и отрезки, ради которых строят график, лежат далеко за её пределами.
Запасы нефти по каждой дате при разных размерах шапки, млн м³
Год
1-й год
m=0,3m = 0{,}3
8,78
m=0,5m = 0{,}5
6,47
m=0,7m = 0{,}7
5,12
Год
2-й год
m=0,3m = 0{,}3
8,75
m=0,5m = 0{,}5
6,47
m=0,7m = 0{,}7
5,13
Год
3-й год
m=0,3m = 0{,}3
8,75
m=0,5m = 0{,}5
6,50
m=0,7m = 0{,}7
5,16
Год
4-й год
m=0,3m = 0{,}3
8,72
m=0,5m = 0{,}5
6,50
m=0,7m = 0{,}7
5,18
Год
5-й год
m=0,3m = 0{,}3
8,68
m=0,5m = 0{,}5
6,49
m=0,7m = 0{,}7
5,19
Год
6-й год
m=0,3m = 0{,}3
8,65
m=0,5m = 0{,}5
6,50
m=0,7m = 0{,}7
5,21

N=F/[(Eo+Efw0)+m(Eg+Efw0)]N = F/\big[(E_o + E_{fw0}) + m(E_g + E_{fw0})\big] по каждой дате. При верном mm запасы от даты не зависят; при заниженном они от даты к дате уменьшаются, при завышенном — растут. Направление дрейфа здесь задаёт нефть урока, у которой Eg/EoE_g/E_o с падением давления уменьшается; у другой нефти оно может быть обратным, но признак тот же — зависимость от даты.

Залежь с водоносным горизонтом

Приток воды WeW_e — не число, а функция времени, и чтобы получить прямую, его выражают через модель водоносного горизонта (Модели притока законтурной воды). Все простые модели имеют вид «параметр горизонта, умноженный на функцию истории давления на контакте»: у малого горизонта функция — само снижение давления Δp\Delta p, у стационарной модели — интеграл снижения давления по времени, у нестационарной — сумма ступеней давления, взвешенных безразмерным притоком (Нестационарный приток воды: метод ван Эвердингена–Херста). Функцию вычисляют по истории давления, а параметр остаётся неизвестным, как и запасы.

Прямая для залежи с водоносным горизонтом
FEt=N+C w(t)Et,Et=Eo+m Eg+Efw\frac{F}{E_t} = N + C\,\frac{w(t)}{E_t},\qquad E_t = E_o + m\,E_g + E_{fw}
где:
  • EtE_tполное расширение на 1 м³ запасов при известном mm (при m=0m = 0 — Eo+EfwE_o + E_{fw}), м³/м³
  • CCпараметр модели водоносного горизонта — наклон прямой, зависят от модели
  • w(t)w(t)функция притока: вычисляемая по истории давления часть модели, We=C w(t)W_e = C\,w(t), зависят от модели
  • N, FN,\ Fгеологические запасы нефти (отрезок прямой) и пластовый объём отбора, м³

Единицы: м³; единицы CC и ww таковы, что их произведение — м³ в пластовых условиях

Справедливо, если
  • тип модели водоносного горизонта и её геометрия выбраны заранее, неизвестен только параметр CC;
  • давление на контакте с водоносным горизонтом принимают равным среднему давлению залежи;
  • размер газовой шапки известен или шапки нет.

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

В других системах единиц (2)
Стационарная модель притока
FEt=N+Ja 1Et∫0t(pi−p) dt\frac{F}{E_t} = N + J_a\,\frac{1}{E_t}\int_0^t (p_i - p)\dd t

JaJ_a — коэффициент продуктивности водоносного горизонта, м³/(сут·МПа) в пластовых условиях; интеграл — МПа·сут; pp — давление на контакте.

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

Малый водоносный горизонт
FEt=N+(cw+cf) W pi−pEt\frac{F}{E_t} = N + (c_w + c_f)\,W\,\frac{p_i - p}{E_t}

WW — поровый объём водоносного горизонта, м³; cw+cfc_w + c_f — сжимаемость воды и пор горизонта, 1/МПа; давление — МПа.

Постоянных нет. Выше давления насыщения EtE_t почти пропорционально Δp\Delta p, абсцисса от даты не меняется, и такой горизонт с помощью графика не выявить.

  • Прямая — модель горизонта согласуется с историей; отрезок — NN, наклон — параметр горизонта.
  • Поздние точки выше прямой, проведённой через ранние — реальный приток растёт быстрее выбранной функции: горизонт больше или проводимее, чем заложено в модель.
  • Поздние точки ниже прямой — приток растёт медленнее: горизонт меньше, давление в нём со временем падает, и он отдаёт воду всё медленнее.
  • S-образная кривая — неверна геометрия, например радиальный горизонт вместо линейного.
  • Беспорядочный разброс — ошибки в данных: давлении, PVT-свойствах, учёте добычи.
  • Точки идут «назад» по оси xx — обычное дело для модели малого горизонта ниже давления насыщения: её абсцисса Δp/Et\Delta p/E_t уменьшается, потому что EtE_t растёт быстрее Δp\Delta p. Выше давления насыщения эта абсцисса почти постоянна, и точки собираются в столбец.
Залежь с водоносным горизонтом: координаты прямой для стационарной модели
Год
1-й год
F/EtF/E_t, млн м³
26,06
∫(pi−p) dt\int(p_i - p)\dd t, МПа·сут
547,9
xx, тыс. МПа·сут
48,77
Год
2-й год
F/EtF/E_t, млн м³
46,47
∫(pi−p) dt\int(p_i - p)\dd t, МПа·сут
2008,9
xx, тыс. МПа·сут
107,10
Год
3-й год
F/EtF/E_t, млн м³
70,20
∫(pi−p) dt\int(p_i - p)\dd t, МПа·сут
4054,3
xx, тыс. МПа·сут
174,83
Год
4-й год
F/EtF/E_t, млн м³
95,44
∫(pi−p) dt\int(p_i - p)\dd t, МПа·сут
6464,9
xx, тыс. МПа·сут
246,95

История — из урока Общее уравнение материального баланса нефтяной залежи: нефтяная часть с начальным давлением 30 МПа, давление всё время выше давления насыщения. x=∫(pi−p) dt/Etx = \int(p_i - p)\dd t/E_t; интеграл — по трапециям между ежегодными замерами давления.

Расчётный примерЗапасы и водоносный горизонт по одной истории
Исходные данные (условные)
История
давление и добыча за четыре года из урока Общее уравнение материального баланса нефтяной залежи (pp = 27,0; 25,0; 23,8; 23,0 МПа при начальном 30 МПа)
Свойства
Boi=1,3243B_{oi} = 1{,}3243; Bw=1,02B_w = 1{,}02 м³/м³; (cwSw+cf)/(1−Sw)=1,2167⋅10−3(c_wS_w + c_f)/(1 - S_w) = 1{,}2167\cdot10^{-3} МПа⁻¹; газовой шапки нет
Запасы
считаются неизвестными; объёмная оценка, 9,06 млн м³, нужна только для сравнения
Решение
  1. Отбор и расширение по годам считаются так же, как в уроке Общее уравнение материального баланса нефтяной залежи; на четвёртый год F=2498,5F = 2498{,}5 тыс. м³, Et=Eo+Efw=0,02618E_t = E_o + E_{fw} = 0{,}02618 м³/м³. Кажущиеся запасы F/EtF/E_t по годам: 26,06, 46,47, 70,20 и 95,44 млн м³. Они растут — в залежь входит вода.
  2. Модель малого горизонта: абсцисса Δp/Et\Delta p/E_t по годам равна 267, 267, 267 и 267 МПа — практически одна и та же, а ординаты различаются почти вчетверо. Прямой через такие точки не провести: приток растёт не только со снижением давления, но и со временем.
  3. Стационарная модель. Интеграл снижения давления — по трапециям, год — 365,25 сут: за первый год 0+32⋅365,25=547,9\tfrac{0 + 3}{2}\cdot365{,}25 = 547{,}9 МПа·сут; к концу второго 547,9+3+52⋅365,25=2008,9547{,}9 + \tfrac{3 + 5}{2}\cdot365{,}25 = 2008{,}9; третьего — 4054,3; четвёртого — 6464,9 МПа·сут.
  4. Абсциссы x=∫(pi−p) dt/Etx = \int(p_i - p)\dd t/E_t: 48,77, 107,10, 174,83 и 246,95 тыс. МПа·сут. Прямая по четырём точкам (метод наименьших квадратов): наклон Ja=350,1J_a = 350{,}1 м³/(сут·МПа), отрезок N=8,98N = 8{,}98 млн м³.
  5. Проверка по двум крайним точкам: Ja=(95,44−26,06)⋅106/[(246,95−48,77)⋅103]=350,1J_a = (95{,}44 - 26{,}06)\cdot10^6/\big[(246{,}95 - 48{,}77)\cdot10^3\big] = 350{,}1 м³/(сут·МПа); N=26,06−350,1⋅48,77⋅10−3=8,99N = 26{,}06 - 350{,}1\cdot48{,}77\cdot10^{-3} = 8{,}99 млн м³.
  6. Приток воды к концу четвёртого года: We=Ja∫(pi−p) dt=350,1⋅6464,9=2263W_e = J_a\int(p_i - p)\dd t = 350{,}1\cdot6464{,}9 = 2263 тыс. м³.
Ответ: Запасы — около 8,98 млн м³, коэффициент продуктивности водоносного горизонта — около 350 м³/(сут·МПа); за четыре года в залежь вошло 2263 тыс. м³ воды.

Что это значит. Объёмная оценка — 9,06 млн м³; расхождение 0,9 % — след округления исходных данных (давлений до 0,1 МПа, объёмного коэффициента нефти до четырёх знаков) и интегрирования по годовым трапециям. В уроке о материальном балансе приток воды нашли, подставив объёмные запасы; здесь те же величины получены из одной истории, без объёмного подсчёта, но ценой предположения о модели водоносного горизонта. Прямая подтверждает стационарную модель лишь на этих четырёх годах: если горизонт конечен, давление в нём со временем снизится, и поздние точки лягут ниже прямой. Поэтому график перестраивают по мере поступления новых замеров. Данные примера условные: история рассчитана с коэффициентом 350 м³/(сут·МПа) и округлена.

Залежь с водоносным горизонтом: точки в координатах стационарной модели и прямая по ним050100150200250020406080100∫(pi − p)dt / Et, тыс. МПа·сутF / Et, млн м³отрезок N = 8,98 млн м³наклон 350 м³/(сут·МПа)

График: по горизонтали абсцисса от 0 до 250 тысяч мегапаскаль-суток, по вертикали F, делённое на Et, от 0 до 100 млн кубометров. Четыре точки — 26 при 49, 46 при 107, 70 при 175 и 95 при 247 — лежат на прямой, которая, продолженная штрихом влево, пересекает ось ординат на уровне 9,0 млн кубометров.

РасчётЧетыре годовые точки залежи с водоносным горизонтом в координатах стационарной модели и прямая, проведённая по ним. Отрезок на оси ординат — запасы нефти, наклон — коэффициент продуктивности водоносного горизонта JaJ_a.

Индексы режимов по истории разработки

Когда модель залежи согласована с историей, отбор можно разделить на вклады источников энергии. Делят на пластовый объём отобранных нефти и газа H=Np[Bo+(Rp−Rs)Bg]H = N_p[B_o + (R_p - R_s)B_g], а из притока воды вычитают добытую воду: место в залежи занимает только та вода, что в ней осталась (индексы режимов).

Индексы режимов и их сумма
NEoH⏟DDI+NmEgH⏟SDI+NEfwH⏟EDI+We−WpBwH⏟WDI=1\underbrace{\frac{NE_o}{H}}_{\text{DDI}} + \underbrace{\frac{NmE_g}{H}}_{\text{SDI}} + \underbrace{\frac{NE_{fw}}{H}}_{\text{EDI}} + \underbrace{\frac{W_e - W_pB_w}{H}}_{\text{WDI}} = 1
где:
  • DDI\text{DDI}индекс расширения нефти (depletion drive index): доля расширения нефти с растворённым в ней газом — режима растворённого газа, а выше давления насыщения — упругого расширения нефти, —
  • SDI\text{SDI}индекс газовой шапки (segregation drive index): доля газонапорного режима, —
  • EDI\text{EDI}индекс упругости воды и пор (expansion drive index): доля расширения связанной воды и уменьшения объёма пор, —
  • WDI\text{WDI}индекс законтурной воды (water drive index): доля водонапорного режима — вошедшей воды за вычетом добытой, —
  • HHпластовый объём отобранных нефти и газа, F−WpBwF - W_pB_w, м³

Единицы: безразмерные доли; все объёмы — в одних единицах

Справедливо, если
  • все величины относятся к одной дате;
  • NN, mm и приток воды найдены по графику или заданы.

Сумма левой части — это [N(Eo+mEg+Efw)+We−WpBw]/H\big[N(E_o + mE_g + E_{fw}) + W_e - W_pB_w\big]/H. Она равна единице ровно тогда, когда точка лежит на прямой; если точка отклоняется от прямой на 0,5 %, сумма отличается от единицы примерно на столько же. Поэтому сумма индексов — ещё и проверка того, насколько хорошо модель описывает каждую дату.

Расчётный примерКак менялся режим двух залежей
Исходные данные (условные)
Залежь с водоносным горизонтом
N=8,98N = 8{,}98 млн м³ и Ja=350,1J_a = 350{,}1 м³/(сут·МПа) — из предыдущего примера; газовой шапки нет
Залежь с газовой шапкой
N=6,5N = 6{,}5 млн м³, m=0,5m = 0{,}5; притока воды нет
Решение
  1. Водоносный горизонт, первый год: H=NpBo=292,8H = N_pB_o = 292{,}8 тыс. м³. Расширение нефти NEo=57,5NE_o = 57{,}5 тыс. м³, вода и поры NEfw=43,4NE_{fw} = 43{,}4 тыс. м³, приток воды Ja⋅547,9=191,8J_a\cdot547{,}9 = 191{,}8 тыс. м³, воды не добывали. Индексы: DDI = 0,196, EDI = 0,148, WDI = 0,655; сумма 1,000.
  2. Четвёртый год: H=2437,3H = 2437{,}3 тыс. м³; приток 2263,4 тыс. м³, добытая вода 60⋅1,02=61,260\cdot1{,}02 = 61{,}2 тыс. м³. Индексы: DDI = 0,055, EDI = 0,042, WDI = 0,904; сумма 1,000.
  3. Газовая шапка, шестой год: H=F=2436,9H = F = 2436{,}9 тыс. м³; нефть NEo=873,8NE_o = 873{,}8, шапка NmEg=1488,6NmE_g = 1488{,}6, вода и поры N(1+m)Efw0=73,0N(1 + m)E_{fw0} = 73{,}0 тыс. м³. Индексы: DDI = 0,359, SDI = 0,611, EDI = 0,030; сумма 0,999.
  4. Газовая шапка, первый год: сумма индексов 1,005 — точка первого года лежит ниже прямой примерно на 0,5 %: при снижении давления всего на 0,53 МПа округление его до 0,01 МПа уже заметно.
Ответ: В залежи с водоносным горизонтом доля законтурной воды выросла за четыре года с 66 до 90 %; в залежи с газовой шапкой около 61 % отбора обеспечивает шапка.

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

Индексы режимов по годам для залежи с водоносным горизонтом и для залежи с газовой шапкойЗалежь с водоносным горизонтом1-й год20 %15 %66 %2-й год11 %81 %3-й год87 %4-й год90 %нефть (DDI)вода и поры (EDI)законтурная вода (WDI)Залежь с газовой шапкой1-й год32 %65 %2-й год32 %64 %3-й год33 %63 %4-й год34 %63 %5-й год35 %62 %6-й год36 %61 %нефть (DDI)газовая шапка (SDI)вода и поры (EDI)

Две группы горизонтальных полос, разделённых на цветные участки. Залежь с водоносным горизонтом, годы с первого по четвёртый: участок законтурной воды растёт с 66 до 90 процентов, участки расширения нефти и воды с порами сокращаются с 20 и 15 до 5 и 4 процентов. Залежь с газовой шапкой, годы с первого по шестой: газовая шапка около 65 процентов в первый год и 61 в шестой, нефть растёт с 32 до 36 процентов, вода и поры — от 4 до 3 процентов.

РасчётИндексы режимов по годам для двух залежей урока: с водоносным горизонтом (запасы и коэффициент горизонта — по прямой из примера) и с газовой шапкой (N=6,5N = 6{,}5 млн м³, m=0,5m = 0{,}5). Длина каждой полосы — сумма индексов; она отличается от единицы не больше чем на половину процента.

Данные, проверки и ограничения

  • Среднее давление решает всё. Расширения — малые разности объёмных коэффициентов при двух давлениях, поэтому ошибка давления переходит в координаты многократно усиленной. Особенно уязвимы ранние точки: координаты делятся на малое расширение. Их не исключают молча, но и не позволяют им определять наклон.
  • PVT-свойства. Объёмные коэффициенты и газосодержание должны отвечать способу разгазирования в пласте и условиям сепарации (Лабораторные PVT-исследования); систематическая ошибка в BoB_o или RsR_s искривляет график так же, как неверная модель.
  • Полный учёт отбора и закачки. Нужен весь добытый газ и вся вода; закачанные вода и газ вычитают из FF. Газ газлифта из добычи исключают.
  • Неоднозначность. Чем больше параметров у модели водоносного горизонта, тем легче «выпрямить» любой график, и тем меньше значит само выпрямление. Предпочитают простейшую модель, согласную с геологией, и пересматривают её с каждым новым замером.
  • Нет времени. График согласует отбор и давление, но не говорит, когда будет достигнуто давление: для прогноза во времени баланс объединяют с расчётом дебитов скважин и притока воды (Прогноз работы залежи с режимом растворённого газа, Модели притока законтурной воды).

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

Метод Хавлены–Одеха
запись уравнения материального баланса как уравнения прямой, по наклону и отрезку которой находят запасы, размер газовой шапки и параметр водоносного горизонта.
FF, EoE_o, EgE_g, EfwE_{fw}
пластовый объём отбора и расширения нефти с её газом, газовой шапки и связанной воды с порами на 1 м³ запасов нефти.
Efw0E_{fw0}
расширение связанной воды и уменьшение объёма пор при m=0m = 0; с газовой шапкой Efw=(1+m)Efw0E_{fw} = (1 + m)E_{fw0}.
Кажущиеся запасы
запасы, рассчитанные по одной дате в предположении, что дополнительных источников энергии нет; их рост во времени указывает на приток воды.
Функция притока w(t)w(t)
вычисляемая по истории давления часть модели водоносного горизонта: We=C w(t)W_e = C\,w(t), где CC — параметр горизонта.
Коэффициент продуктивности водоносного горизонта JaJ_a
параметр стационарной модели: приток воды за сутки на 1 МПа снижения давления на контакте.
Индексы режимов (DDI, SDI, EDI, WDI)
доли пластового объёма отобранных нефти и газа, обеспеченные расширением нефти, газовой шапкой, водой с порами и законтурной водой.

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

Неверно

«Точки легли на прямую — значит, модель залежи верна».

На самом деле

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

Неверно

«По графику F/Eo от Eg/Eo можно найти и запасы, и размер газовой шапки».

На самом деле

Можно только при длинной истории и очень точном давлении. Абсциссы точек меняются мало, отрезок получается долгой экстраполяцией, и округление давления до 0,1 МПа в примере урока превратило 6,5 млн м³ в 1,6, а без самой чувствительной, первой точки — в 7,6. Размер шапки берут по геологическим данным, а баланс служит проверкой.

Неверно

«Точки изогнулись — надо подобрать другой коэффициент водоносного горизонта».

На самом деле

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

Неверно

«Первые замеры самые надёжные: залежь ещё не тронута».

На самом деле

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

Неверно

«Сумма индексов режимов не равна единице — значит, ошибка в арифметике».

На самом деле

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

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

Та же прямая для газовой залежи

Для газовой залежи с водоносным горизонтом отбор F=GpBg+WpBwF = G_pB_g + W_pB_w равен расширению газа G(Bg−Bgi)G(B_g - B_{gi}) и притоку воды; расширением связанной воды и пор обычно пренебрегают (Материальный баланс газовой залежи). Деление на расширение кубометра газа даёт прямую того же вида.

Прямая для газовой залежи со стационарным водоносным горизонтом
GpBg+WpBwBg−Bgi=G+Ja 1Bg−Bgi∫0t(pi−p) dt\frac{G_pB_g + W_pB_w}{B_g - B_{gi}} = G + J_a\,\frac{1}{B_g - B_{gi}}\int_0^t (p_i - p)\dd t
где:
  • G, GpG,\ G_pгеологические запасы газа (отрезок прямой) и накопленная добыча газа в стандартных условиях, м³
  • Bg, BgiB_g,\ B_{gi}объёмный коэффициент газа при текущем и начальном давлении, м³/м³
  • Wp, BwW_p,\ B_wдобытая вода и её объёмный коэффициент, м³, м³/м³
  • JaJ_aкоэффициент продуктивности водоносного горизонта — наклон прямой, м³/(сут·МПа)
  • pi, pp_i,\ pначальное и текущее давление на контакте (среднее пластовое), абсолютные, МПа

Единицы: м³ в стандартных условиях; интеграл — МПа·сут

Справедливо, если
  • стационарная модель водоносного горизонта;
  • расширением связанной воды и пор пренебрегают; газ в пласте однофазен.

Левая часть — кажущиеся запасы газа GaG_a урока о газовом балансе. Там они росли от замера к замеру, и это было лишь признаком притока воды; здесь их рост объяснён количественно, и запасы находятся вместе с параметром горизонта. Расширение кубометра газа записано явно, Bg−BgiB_g - B_{gi}: символ EgE_g в нефтяном балансе означает другую величину.

Расчётный примерГазовая залежь: запасы и водоносный горизонт вместе
Исходные данные (условные)
История
залежь с водоносным горизонтом из урока Материальный баланс газовой залежи: пять ежегодных замеров давления, отбор 1,0 млн м³/сут, воду не добывают
Начальные условия
pi=20p_i = 20 МПа, Bgi=0,005302B_{gi} = 0{,}005302 м³/м³; запасы считаются неизвестными
Решение
  1. Для каждого года: y=GpBg/(Bg−Bgi)y = G_pB_g/(B_g - B_{gi}) — кажущиеся запасы. Например, на пятый год p=14,62p = 14{,}62 МПа, Bg=0,007210B_g = 0{,}007210 м³/м³, y=6,90y = 6{,}90 млрд м³; в первый год — 5,38 млрд м³.
  2. Интеграл снижения давления по трапециям: к концу пятого года 5464 МПа·сут; абсциссы x=∫(pi−p) dt/(Bg−Bgi)x = \int(p_i - p)\dd t/(B_g - B_{gi}) — от 0,69 до 2,86 млн МПа·сут.
  3. Прямая по пяти точкам: отрезок G=4,91G = 4{,}91 млрд м³, наклон Ja=695,5J_a = 695{,}5 м³/(сут·МПа).
Ответ: Запасы газа — около 4,91 млрд м³ при объёмной оценке 4,89, коэффициент продуктивности водоносного горизонта — около 696 м³/(сут·МПа).

Что это значит. В уроке о газовом балансе приток воды нашли по объёмным запасам; здесь запасы и водоносный горизонт определены вместе. История там рассчитана с коэффициентом 700 м³/(сут·МПа) по месячным шагам, поэтому годовые трапеции и округление давлений до 0,01 МПа дают небольшое расхождение. Данные примера условные.

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

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

  1. 1

    Кажущиеся запасы залежи без газовой шапки, рассчитанные по четырём годам, — 9,1; 9,0; 9,1 и 9,0 млн м³; давление всё время выше давления насыщения, объёмный подсчёт — 6,0 млн м³. Что это может означать?

    1. А.Ошибка в накопленном газовом факторе
    2. Б.Залежь замкнута, объёмный подсчёт верен
    3. В.В залежь входит вода из большого водоносного горизонта
    4. Г.Есть дополнительный источник энергии, меняющийся с давлением так же, как расширение нефти, — например, небольшой водоносный горизонт, — либо занижен объёмный подсчёт
    Показать ответ

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

    Постоянство кажущихся запасов говорит, что дополнительный источник энергии, если он есть, пропорционален Δp\Delta p, как и расширение недонасыщенной нефти. Небольшой водоносный горизонт выше давления насыщения ведёт себя именно так и поднимает уровень линии, не меняя её формы. Большой горизонт дал бы рост кажущихся запасов со временем. Выбрать между горизонтом и ошибкой объёмного подсчёта поможет геология законтурной области.

  2. 2

    Две точки графика F/EoF/E_o — Eg/EoE_g/E_o: x=4,0x = 4{,}0, y=15,0y = 15{,}0 млн м³ и x=3,0x = 3{,}0, y=12,5y = 12{,}5 млн м³. Найдите запасы нефти и размер газовой шапки. Как изменится результат, если ордината первой точки на 1 % больше?

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

    Наклон mN=(15,0−12,5)/(4,0−3,0)=2,5mN = (15{,}0 - 12{,}5)/(4{,}0 - 3{,}0) = 2{,}5 млн м³, отрезок N=12,5−2,5⋅3,0=5,0N = 12{,}5 - 2{,}5\cdot3{,}0 = 5{,}0 млн м³, m=2,5/5,0=0,5m = 2{,}5/5{,}0 = 0{,}5. При y1=15,15y_1 = 15{,}15: наклон 2,65, N=12,5−2,65⋅3,0=4,55N = 12{,}5 - 2{,}65\cdot3{,}0 = 4{,}55 млн м³ (на 9 % меньше), m=0,58m = 0{,}58. Ошибка в 1 % выросла в девять раз: при продолжении прямой до оси ординат абсолютная ошибка утроилась — с 0,15 до 0,45 млн м³, потому что отрезок лежит втрое дальше от точек, чем они друг от друга, — а сам отрезок втрое меньше ординаты точки.

  3. 3

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

    1. А.В залежи появилась газовая шапка
    2. Б.Реальный приток растёт медленнее, чем интеграл снижения давления: горизонт конечен, и давление в нём падает
    3. В.Нужно увеличить коэффициент JaJ_a
    4. Г.Водоносный горизонт больше, чем предполагалось
    Показать ответ

    Правильный ответ: Б. Реальный приток растёт медленнее, чем интеграл снижения давления: горизонт конечен, и давление в нём падает

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

  4. 4

    На некоторую дату: H=2000H = 2000 тыс. м³, добытая вода WpBw=100W_pB_w = 100 тыс. м³; NEo=300NE_o = 300, NmEg=500NmE_g = 500, NEfw=60NE_{fw} = 60, We=1240W_e = 1240 тыс. м³. Найдите индексы режимов и их сумму.

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

    DDI = 300/2000 = 0,15; SDI = 500/2000 = 0,25; EDI = 60/2000 = 0,03; WDI = (1240 − 100)/2000 = 0,57. Сумма — 1,00: отбор F=H+WpBw=2100F = H + W_pB_w = 2100 тыс. м³ точно равен сумме источников, 300+500+60+1240=2100300 + 500 + 60 + 1240 = 2100 тыс. м³, то есть точка лежит на прямой. Больше половины отбора дала законтурная вода, четверть — газовая шапка.

  5. 5

    Почему запасы нефти и размер газовой шапки трудно определить вместе по одному балансу?

    1. А.Расширения нефти и газа шапки растут с падением давления почти пропорционально, поэтому точки на графике собираются в узкую полосу и отрезок находят долгой экстраполяцией
    2. Б.Газ шапки не входит в уравнение баланса
    3. В.Размер шапки меняется во времени
    4. Г.Баланс неприменим ниже давления насыщения
    Показать ответ

    Правильный ответ: А. Расширения нефти и газа шапки растут с падением давления почти пропорционально, поэтому точки на графике собираются в узкую полосу и отрезок находят долгой экстраполяцией

    Отношение Eg/EoE_g/E_o за историю меняется мало: в примере урока — от 4,1 до 3,4 (абсцисса с учётом связанной воды и пор — от 3,8 до 3,3). Отрезок на оси ординат лежит далеко от точек, и ошибки давления в несколько сотых мегапаскаля сдвигают его на десятки процентов. Надёжно определяется только комбинация N(1+mxˉ)N(1 + m\bar x); поэтому mm берут по геологическим данным.

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