Научная статья на тему 'Метод материального баланса как первичный инструмент оценки показателей разработки участка месторождения при заводнении'

Метод материального баланса как первичный инструмент оценки показателей разработки участка месторождения при заводнении Текст научной статьи по специальности «Науки о Земле и смежные экологические науки»

CC BY
3086
308
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
УРАВНЕНИЕ МАТЕРИАЛЬНОГО БАЛАНСА / ОТКРЫТЫЙ УГОЛ ПРИТОКА / ГЕОМЕТРИЧЕСКИЙ КОЭФФИЦИЕНТ УЧАСТИЯ / ПРИТОК ИЗ ЗАКОНТУРНОЙ ОБЛАСТИ / MATERIAL BALANCE EQUATION / ANGLE OPEN-TO-FLOW / GEOMETRICAL WELL ALLOCATION FACTOR / WATER INFLUX

Аннотация научной статьи по наукам о Земле и смежным экологическим наукам, автор научной работы — Абидов Денис Геннадьевич, Камартдинов Марат Рамисович

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по наукам о Земле и смежным экологическим наукам , автор научной работы — Абидов Денис Геннадьевич, Камартдинов Марат Рамисович

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

The article considers the problem of applying the material balance when estimating the efficiency of field area flooding. The authors propose the solution of the problem of distributing the produced and pumped liquid volumes of total edge wells between the adjacent areas (the problem of well allocation estimation). The method is based on transition from constant geometrical well allocation factors which are valid only at symmetric pattern of pressure distribution to their alternating analogues. It could reflect the real situation when the pressure distribution pattern changes. This effects, in its turn, on well flow rates distribution among the areas.

Текст научной работы на тему «Метод материального баланса как первичный инструмент оценки показателей разработки участка месторождения при заводнении»

Геология нефти и газа

УДК 622.276.43

МЕТОД МАТЕРИАЛЬНОГО БАЛАНСА КАК ПЕРВИЧНЫЙ ИНСТРУМЕНТ ОЦЕНКИ ПОКАЗАТЕЛЕЙ РАЗРАБОТКИ УЧАСТКА МЕСТОРОЖДЕНИЯ ПРИ ЗАВОДНЕНИИ

Д.Г. Абидов*, М.Р. Камартдинов

Томский политехнический университет E-mail: kamartdinovmr@hw.tpu.ru *Компания «Сахалин Энерджи Инвестмент Компани Лтд.», г. Южно-Сахалинск E-mail: Denis.Abidov@sakhalinenergy.ru

Рассмотрена проблема применения материального баланса взадаче оценки эффективности заводнения участка месторождения. Предложено решение проблемы распределения объемов добытой изакачанной жидкости отобщих граничных скважин между соседними участками (проблема оценки участия скважин). Метод основывается на переходе от постоянных геометрических коэффициентов участия скважин, корректных лишь в условиях симметричной картины распределения давления, к ихпе-ременным аналогам, что отражало бы реальную ситуацию, когда картина распределения давления меняется, что в свою очередь влияет на распределение дебитов скважин между участками.

Ключевые слова:

Уравнение материального баланса, открытый угол притока, геометрический коэффициент участия, приток из законтурной области.

Key words:

Material balance equation, angle open-to-flow, geometrical wellallocation factor, water influx.

Уравнение материального баланса (МБ) - универсальный инструмент, позволяющий оценить запасы углеводородов, объемы воды, поступающей из законтурной области, и режимы разработки месторождения. Очевидное преимущество использования уравнения МБ по сравнению с процессом построения и адаптации гидродинамической модели месторождения - использование минимального набора исходной информации [1].

В классическом представлении теории МБ пласт рассматривается как изолированный по-ровый объем, который на любой момент времени можно охарактеризовать величиной среднего пластового давления. Применение уравнения МБ для участка месторождения рассматривалось несколькими авторами [2, 3]. Но рассматриваемый в [2, 3] подход основан на использовании постоянных геометрических коэффициентов участия скважин, корректных лишь в условиях симметричной картины распределения давления.

Целью данной работы является расширение традиционной концепции применения уравнения МБ для анализа эффективности процесса заводне-

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

В качестве апробации предложенной методики в данной статье рассматривается одно из нефтяных месторождений Томской области, на котором сформирована система поддержания пластового давления - пластовое давление поддерживается выше давления насыщения. Следовательно, уравнение МБ принимает следующий вид:

NPB0 + WpBw -WiBw =

N (Bo - Boi) + N (cf + cAc) AP + W, (1)

1 Swc

где Np - накопленная добыча нефти, м3; N - начальные геологические запасы нефти, м3; Wp - накопленная добыча воды, м3; W - накопленный объем закачанной воды, м3; Вы, Во - объемный коэффициент нефти при начальном (26 МПа) и при

текущем пластовом давлении, м3/м3; Вк - объемный коэффициент воды, м3/м3; Зж - насыщенность связанной воды, д.е.; сг - сжимаемость поро-вого пространства, атм-1; е„ - сжимаемость воды, атм-1; АР=Р-Р - снижение давления, разница между начальным пластовым давлением и величиной среднего пластового давления на текущий момент, атм; Же - накопленный приток воды из законтурной области, м3.

Анализируемое месторождение имеет сложное геологическое строение. Нефтенасыщенными продуктивными интервалами являются пласты юрской группы Ю11Л, Ю11, Ю10. Каждый из пластов имеет определенную протяженность и выдержанность по площади месторождения. Можно выделить три основных типа разреза, представленных: 1 - относительно тонким по мощности песчаным пластом; 2 - двумя пластами, разделенными глинистой перемычкой; 3 - сравнительно мощным песчаным пластом, в котором присутствуют тонкие непроницаемые пропластки. На месторождении сформирована трехрядная схема заводнения (с тремя рядами добывающих скважин и одним рядом нагнетательных), рис. 1. Месторождение разделено рядами нагнетательных скважин на 17 блоков.

Так как результат применения уравнения МБ сильно зависит от физико-химических свойств пластовых флюидов [4], то для анализа был выбран блок № 15 в той части месторождения, где скважины вскрывают только один пласт Ю10, рис. 1. Это позволит избежать неопределенностей, связанных с оценкой физико-химических свойств пластовых флюидов, а также с оценкой средних значений фильтрационно-емкостных свойств пластов.

При анализе уравнения МБ основной исходной информацией являются данные по накопленной добыче жидкости и закачке воды. При попытке применить уравнение МБ на выбранном элементе разработки важно помнить, что не весь объем добываемого или закачиваемого скважинами флюида относится именно к выбранному блоку. Часть флюидов может поступать от соседних блоков. Таким образом, для использования уравнения МБ необходимо прежде всего распределить добычу и закачку по анализируемым элементам разработки.

Рис. 1. Карта проводимости условного нефтяного месторождения

Корректное распределение добываемых и закачиваемых объемов жидкости сильно зависит от распределения давления вокруг скважин и фактически возможно только после построения полноценной гидродинамической модели месторождения [3]. В данной статье предлагается новый подход к оценке распределения добываемых и закачиваемых объемов жидкости, в котором правильность полученного распределения оценивается с помощью уравнения МБ.

Один из простейших подходов состоит в том, что объемы определяются пропорционально углам притоков в каждый элемент, имеющий общую скважину [5]. Два простейших примера иллюстрируют этот метод (рис. 2). Согласно рис. 2, а, угловые скважины нагнетают в выбранный элемент (ограниченный красной линией) только 120/360=1/3 от общего объема закачки, тогда как для ситуации, представленной на рис. 2, б, доля закачиваемых в элемент объемов составляет 90/360=1/4. Коэффициент, рассчитываемый на основе угла притока, называется геометрическим коэффициентом участия и обозначается как WAF (well allocation factor).

\ \ / / g ^ А к • 4 «--1 1 • : 1

t>~* ■ Ь • i *\90° ^ • i 1 А

» » 1 щ m ' ' / \ ▲ -A i--4 a i • 1 Г Т 1 i • À

Рис. 2. Оценка распределения флюидов по блокам, основанная на углах притока. Схема размещения скважин: а) семиточечная; б) девятиточечная

Данный подход, использующий угол притока для распределения объемов жидкости по блокам, подходит для геологически однородных пластов, с симметричными схемами расстановки скважин, работающих с одинаковыми дебитами при одинаковых забойных давлениях. Однако на практике такие условия практически никогда не выполняются, что приводит к нарушению симметрии и варьированию во времени величины ЩА!, что можно продемонстрировать на примере модели линий тока [6], рис. 3.

Таким образом, можно сделать вывод, что используемое определение величины Ш! неподхо-дит для решения поставленной задачи и необходимо использовать аналогичную переменную, являющуюся функцией времени.

Поскольку в ур. (1) величины Лр, Щ Щявляют-ся накопленными объемами, то они не могут быть использованы напрямую для анализа выбранного элемента разработки. Для того чтобы учесть распределение добываемых и закачиваемых объемов по блокам, необходимо, чтобы объемы добычи и закачки были правильно распределены между элементами разработки:

N * В + Ж * В =

р о р п

= N (Во - Вт) + N

В

~(сг + €^с) + ЖВ , (2)

где Л/, Щ/, Щ* - соответствующие объемы с учетом распределения добычи и закачки по блокам.

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

ап,_ = •

360'

(3)

ак (*) = а0к + ак (()•

(4)

В ур. (4) а0к - величина известная, согласно определению угла притока (рис. 1 и 2). Величину а0к можно использовать как первое приближение при оценке ак(1), а величина ак(1) позволяет оценить более точное значение ак(1) и его зависимость от времени. Следовательно, распределенные объемы добычи нефти будут рассчитываться согласно следующему выражению:

Чк({) = ак 0 И 0 )•

Добыча и закачка воды рассчитывается аналогичным образом.

Таким образом, если анализируемый блок содержит N скважин (в нашем случае 15 блок содержит 94 скважины), то необходимо определить N неизвестных ак(1) на каждый временной шаг. Но при этом мы можем использовать только одно уравнение - уравнение МБ. Очевидно, что необходимо сократить количество переменных, заменяя N коэффициентов ак(1) одним общим коэффициентом ав(1), который можно представить как среднее арифметическое всех ак(1):

а в (() =

ВОРМ *(;) ВОРМ (Г)

Хак 0) • Чк 0)

к =1______________

N -

X Ук 0)А^

(5)

Для учета изменения ЩА! со временем введем дополнительную переменную, зависящую от времени, таким образом, чтобы динамический ЩА! принял вид:

где ВОРМ (?) = Х Ук ({)^ - месячная (если вре-

к =1

менной шаг анализа составляет 1 месяц) накопленная добыча нефти по анализируемому блоку без учета распределения объемов флюидов, а

N

ВОРМ*(1) = ^ак({)Чк({)^ - месячная накоплен-

к=1

ная добыча нефти по анализируемому блоку с учетом распределения объема. Аналогичным образом определяются месячная накопленная добыча и за-

к=1

качка воды ВЩРМ(1), ВЩ1М(1) и соответствующие коэффициенты распределения объемов добычи и закачки воды по блокам.

Совокупный коэффициент распределения объемов добычи и закачки для всего блока определяется аналогично коэффициенту ак(1):

ав (0 = ав0 +(ав (0> (6)

где аВ0(1) - геометрический коэффициент участия для всех скважин блока, введенный по аналогии с коэффициентом а0к (3), который можно оценить из выражения (5), приравнивая в ур. (4) ак(1) нулю.

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

N *

Рв ^)=-Е-

N

£ав (г) ВОРМ (г) _0______________

г '

X ВОРМ (г)

(7)

Совокупный коэффициент распределения накопленных объемов для всего блока определяется как:

Рв (г) = Рв 0 + Рв (г )•

Таким образом, количество неизвестных было сокращено до одной неизвестной величины аВ(1), определив которую, с помощью ур. (7) можно определить и Рв(1).

На каждый временной шаг в ур. (2) используются следующие распределенные объемы добычи и закачки (индексы определяют «ро» - добыча нефти; «рт - добыча воды; «т»> - нагнетание воды):

N..

Ж* = ЖРп,

К = Жр Ррп •

В общем случае, при анализе выбранного блока месторождения, величина накопленных объемов притока воды из законтурной области Щ также является неизвестной. Таким образом, в ур. (2) имеем 4 неизвестных величины на каждый временной шаг авр„(1), авр/1), ат(1), Щ. Допуская, что до и после распределения по блокам объемов добычи жидкости обводненность в скважинах не изменяется, аВрт(1)=аВр,0(1). Таким образом, число неизвестных в ур. (2) сократилось до трех.

Величина Щ традиционно рассчитывается с помощью уравнения МБ, но для этого необходимо, чтобы были известны накопленные объемы добычи и закачки жидкости. Аналогичным подходом воспользуемся и в нашем случае. Предположим, что в течение некоторого периода Хоценка распределения объемов добычи и закачки с помощью геометрического НАР дает достаточную точность. Следовательно, можно считать, что величины Лр*, Щр*, Щ* известны. Используя различные модели расчета притока из законтурной области Щ, добиваемся равенства левой и правой части ур. (2) и по-

лучаем основные параметры и свойства законтурной области. Основная проблема в данном подходе - достоверное выделение периода Х

Начальный период

■ IV, (Аналитический)

а

Время, год

- IV, (МБ)

■ IV' (Аналитический)

Рис. 4.

б

Сравнение величины притока из законтурной области. Результаты расчета: а) исходные; б) после подбора коэффициентов аВро(1) и ав,М

Из анализа информации об истории разработки рассматриваемого месторождения известно, что ввод скважин в работу в выбранном блоке производился постепенно в направлении с северо-запада на юго-восток. Было сделано предположение, что наиболее достоверным периодом для определения параметров законтурной области являются первые два года разработки выбранного блока; во внимание были приняты следующие факторы:

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

2. Если в этот период и существуют неучтенные объемы перетоков жидкости в пределах блока, то данные объемы должны в большей степени определяться притоком из законтурной области, так как количество пробуренных скважин сравнительно невелико.

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

Предполагаем, что начальный период добычи дает меньшую накопленную ошибку при ра-

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

3.

4.

0

спределении накопленных объемов добычи и закачки для выбранного блока.

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

Величина Щ рассчитывается по методу МБ (2) и с помощью аналитической модели притока из законтурной области. Параметры аналитической модели подбираются таким образом, чтобы результаты расчетов двумя методами согласовались в течение первых двух лет (предположение, что распределение объемов добычи и закачки вычисляется на основе геометрического коэффициента участия, который постоянен во времени в этот период). Как и предполагалось, на более поздних этапах различия увеличиваются (рис. 4, а) - это связанно с тем, что коэффициент участия скважин меняется со временем.

Таким образом, получив параметры аналитической модели для расчета Щ, количество неизвестных в ур. (2) сократилось до двух: аВро(1) и аВм(1). При этом предполагаем, что параметры для расчета притока из законтурной области являются истинными и не будут меняться со временем. Задача сводится к подбору таких значений аВря(1) и аВм(1), чтобы минимизировать различия в расчете Щ по методу МБ и с помощью аналитической модели (рис. 4, а). Результат показан на рис. 4, б, разница в расчете Щ была устранена за счет распределения объемов добычи и закачки флюидов по блокам с помощью коэффициентов аВро(1) и аВм(1).

Рис. 5. Коэффициенты распределения накопленных объемов добычи нефти и закачки воды

Согласно полученным результатам подбора коэффициентов ав,ро(0 и ав,ы(1), а следовательно и значений Д/І) и Ды(1) (рис. 5), выделим временные периоды:

1. От 0 до 2 лет: ДД^Д^О, Д,(0<Д,о(0. Из анализируемого блока добывается больше нефти, чем это предполагалось на основе геометрического коэффициента участия, а воды закачивается меньше. Данный факт можно объяснить следующими факторами: соседний блок № 14 к

тому времени уже активно разрабатывается, и пластовое давление в нем ниже, чем давление в анализируемом блоке № 15. Таким образом, добыча нефти скважинами идет в основном за счет притока нефти в пределах блока с более высоким пластовым давлением (блок № 15), а закачка, наоборот, направляется в зону пониженного давления - в соседний блок № 14.

2. От 2 до 10 лет: Д^КД^О, РМ>Рм(1). Из анализируемого блока добывается меньше нефти, чем это предполагалось на основе геометрического коэффициента участия, а воды закачивается больше. В этот период практически все скважины блока уже введены в эксплуатацию. Нефтенасыщенность в пределах блока существенно снизилась, что предполагает увеличение доли добываемой нефти, поступающей в добывающие скважины из-за пределов анализируемого блока. Снижение пластового давления в блоке является причиной того, что нагнетаемая вода остается в пределах анализируемого участка.

3. От 10 лет: РДО^ДоЛО, ДМ^Рм(). Добывающие скважины находятся внутри блока и не граничат с соседними участками, т. е. добыча определяется в основном притоком нефти в пределах блока. Коэффициент распределения объемов добычи нефти приближается к постоянному значению, соответствующему величине геометрического коэффициента участия (который близок к 1). Распределение закачиваемых в пласт объемов воды соответствует геометрическому коэффициенту участия, так как смежные блоки 14-16 характеризуются сравнительно одинаковыми пластовыми давлениями и закачка в блоки происходит пропорционально углу притока.

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

Выводы

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

СПИСОК ЛИТЕРАТУРЫ

1. Дейк Л.П. Практический инжиниринг резервуаров. - М.: Институт компьютерных исследований, 2008. - 668 с.

2. Уолкотт Д.С. Разработка и управление месторождениями при заводнении. - М.: Юкос, 2001. - 144 с.

3. Ehlig-Economides С.А. Applications for Multiphase Compartmentalized Material Balance // University of Tulsa: Centennial Petroleum Engineering Symposium. - Tulsa, 1994. - V. 1. - P. 429-441.

4. Baker R.O. Reservoir management for waterfloods. P. II // Journal of Canadian Petroleum Technology. - 1998. - V. 37. - № 1. -P. 12-17.

5. Chapman L.R., Thomson R.R. Waterflood Surveillance in the Ku-paruk River Unit with Computerized Pattern Analysis // Journal of Petroleum Technology. - 1989. - V. 41. - № 3. - P. 277-282.

6. Chugh S.H., Baker R.O. Revisiting reservoir flood-surveillance methods using streamlines // SPE Reservoir Evaluation & Engineering. - 2008. - V. 11. - № 2. - P. 387-394.

Поступила 17.04.2012 г.

УДК 551.762(571.1)

СХЕМА ИНДЕКСАЦИИ ПРОДУКТИВНЫХ ПЕСЧАНЫХ ПЛАСТОВ ЮРЫ ЗАПАДНОЙ СИБИРИ НА ОСНОВЕ ИНВЕРСИОННОЙ МОДЕЛИ ЦИКЛОГЕНЕЗА

А.Л. Бейзель

Институт нефтегазовой геологии и геофизики им. А.А. Трофимука СО РАН, г. Новосибирск E-mail: beiselal@ipgg.nsc.ru

Разработана новая схема индексации песчаных продуктивных пластов, основанная на раздельном подходе к континентальным и морским разрезам. В аллювиальных циклах песчаные пласты базальные, а в морских - кровельные, регрессивные. Они не могут быть синхронными друг другу. В предлагаемой схеме вместо единой серии индексов вида Юп представлены две колонки индексов: морские пласты обозначены буквой «м», а континентальные - «к», например: Юм1 иЮк1. Для пластов неопределенной фациальной принадлежности идля циклов в целом остается индекс вида Ю. Пласты типа Ю20 рассматриваются как непосредственное продолжение континентальных базальных пластов главной фазы активизации речного стока. Самый нижний пласт васюганского горизонта в морских фациях будет иметь индекс Юм20, а синхронный ему континентальный - Юк20.

Ключевые слова:

Юра, Западная Сибирь, индексация пластов.

Key words:

Jurassic, Western Siberia, productivity layers indexing.

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

Индексация представлена в региональных стратиграфических схемах, где пласты или их группы отмечены в характеристиках свит. Этим как бы подчеркивается их неопределенное стратиграфическое положение внутри литостратонов. Для юры Западной Сибири схемы индексации в наиболее общем виде представлены в работах [1, 2]. Сразу следует отметить, что они сделаны независимо от фациальной принадлежности отложений. Морские песчаные пласты оказались увязанными с континентальными и получили одинаковые ин-

дексы. Между тем, в разнофациальных отложениях не принято выделять даже единые свиты.

Многими исследователями было отмечено, что в терригенных отложениях в зависимости от обстановок осадконакопления (континентальной или морской) имеет место смена гранулометрического состава пород разной направленности. Для континентальных обстановок характерно уменьшение размерности обломочных зерен вверх по разрезу, в то время как для отложений морского генезиса преобладающей является обратная последовательность распределения обломочных зерен [3-8]. В Западной Сибири хорошо известна приуроченность проциклитов к континентальным отложениям нижней и средней юры, а рециклитов и асимметричных прорециклитов - к морским осадкам келловея и верхней юры [9-11]. Морские песчаные пласты являются регрессивными и формируются в конце цикла, а континентальные базальные пласты - в начале циклов. Они не могут быть синхронными друг другу. Если предположить, что циклогенез в разных фациальных областях идет независимо, то еще менее вероятно, что песчаные пласты в них будут совпадать между собой.

i Надоели баннеры? Вы всегда можете отключить рекламу.