пластов по методике ВНИИ (методика Ю.П. Борисова).
Ю.П. Борисовым предложен метод расчета показателей разработки нефтяных залежей с учетом неоднородности пластов по проницаемости. Метод Борисова является основой для применяемых в настоящие время методов расчета процесса разработки нефтяных залежей при заводнении.
В качестве расчетной модели принимается пласт, состоящей из набора параллельно работающих трубок тока одинакового поперечного сечения. Трубки тока имеют разную проницаемость и вероятностно распределены в объеме пласта.
В первом приближении принято, что число трубок тока определяется эмпирической кривой распределения проницаемости, построенной по данным кернового или геофизического материала.
Предполагается, что расход жидкости через трубку тока в каждый момент времени пропорционален ее проницаемости.
Ряды скважин рассматриваются как эквивалентные галереи с дополнительным внутренним фильтроционным сопротивлением.
Очень важный фактор, что в действительности не происходит полного, поршневого вытеснения нефти водой и за фронтом вытеснения остается водонефтяная смесь, нефтенасыщенность которой уменьшается по мере прокачки жидкости, учитывается преобразованием спектра распределения трубок тока.
Кривая изменения насыщенности при вытеснении нефти водой:
 |
где: z – насыщенность породы подвижной нефтью в зоне водонефтяной смеси;
z ф – на фронте вытеснения;
S с.в. – количество связанной воды или первоначальная водонасыщенность;
S о.н. – остаточная нефтенасыщенность или количество нефти, которая остается
породе после бесконечно долгой промывки;
S н.ф. – общая нефтенасыщенность на фронте вытеснения;
V ф – объем пласта, до фронта продвижения ВНК.
Непоршневая часть эпюры вытеснения рассматривается как объем пласта, в котором трубки тока обладают дополнительной неоднородностью по проницаемости.
Эта проницаемость изменяется по закону, соответствующему распределению насыщенности на данном участке:
(1)
где: z – насыщенность подвижной нефтью;
m , V – пористость и объем пласта;
– количество вторгшейся в пласт жидкости.
Уравнение (1) получено Борисовым на основании обработки кривых фазовых проницаемостей Эфроса и уравнения Баклея – Леверетта и справедливо для соотношения 1≤ 
≤10
Насыщенность на фронте вытеснения z ф определяется из уравнения:
(2)
Величину z ф находят методом итераций (последовательных приближений).
Схема расчета процесса вытеснения из преобразованного пласта строится путем определения характеристик вытеснения по отдельным трубкам тока, где вытеснение считается поршневым. Затем суммируются результаты по всем трубкам с различной проницаемостью.
Дальнейшие расчеты технологических показателей разработки выполняются или при заданных перепадах давления, или при заданных дебитах жидкости для различных систем заводнения.
Для однорядной системы заводнения при условии равнодебитности эксплуатируемых и нагнетательных скважин, дебит одной скважины или перепад давления определяется из уравнения:
(3)
где rc , r н – радиус добывающей и нагнетательной скважины;
r ф– радиус фронта вытеснения;
d – расстояние от нагнетательной скважины до ближайшей эксплуатационной
скважины.
Время для соответствующих положений определяется по формуле:
(4)
Для девятиточечной площадной системы заводнения дебит одной нагнетательной скважины или трех эксплуатационных или перепад давления определяется из уравнения:
(5)
Время для различных положений при Δ P = const ( 
) определяется по формуле:
(6)
где 
– средний коэффициент использования объема пор в пределах водонефтяной
зоны.
При проектировании разработки неоднородных прерывистых пластов технологические показатели предлагается рассчитывать двумя способами: по схеме прерывистого пласта и по схеме непрерывистого пласта с использованием дополнительных коэффициентов, учитывающих прерывистость и линзовидность.
По схеме прерывистых пластов объем схематизируется серией объемов линз и полулинз различной длины. Данные о протяженности линз и полулинз по отношению к направлению фильтрации жидкости получают в результате обработки зональных карт распространения прослоев и составление при этом ранжированного ряда длин и объемов линз и полулинз. Из расчета исключаются все объемы линз.
Расчетной схемой является набор пластов различной длины, при этом начальное положение ВНК для всех пластов принимается одинаковым.
Расчеты технологических показателей разработки неоднородных пластов по методу УФНИИ (метод М.М. Саттарова).
М.М. Саттаровым предложен метод расчетов показателей разработки с учетом неоднородности по проницаемости.
Метод предполагает, что непрерывный неоднородный пласт моделируется серией прослоев, границы которых сориентированы параллельно линиям тока. Прослои различной проницаемости распределены по мощности вероятностно, а плотность f ( k ) и функция распределения проницаемости F ( k ) имеет вид видоизмененного распределения Максвелла:
(1)
(2)
где k – проницаемость,
a , 
– параметры распределения.
Вытеснения нефти принимается поршневым, но учитывается снижение фазовой проницаемости для воды в промытой зоне.
Принимается течение жидкости к “проницаемым” эксплуатационным галереям с дополнительным внутренним фильтроционным сопротивлением, а скорость движения жидкости – пропорциональной проницаемости пропластков.
Схема расчета процесса обводнения строится аналогично схеме расчетов по методу Ю.П. Борисова – определяются характеристики вытеснения по отдельным прослоям, которые затем суммируются по всем прослоям различной проницаемости.
Доля нефти в потоке жидкости находится по следующей формуле:
(3)
где: 
– вязкость нефти и воды,
– относительная проницаемость для воды в заводненной зоне (принимается
равной 0,5÷0,3 по кривым фазовой проницаемости Березина В.М.),
– средняя проницаемость зоны, занятой вытесняемой водой к моменту tm,
– средняя проницаемость зоны, занятой нефтью к моменту tm .
С учетом функции распределения проницаемости (2) 
и 
находятся из соотношений:
(4)
(5)
Коэффициент охвата пласта заводнением к моменту tm прорыва воды по пропласткам с проницаемостью km представляется как:
(6)
где: F ( km ) – значение интегральной функции распределения проницаемости в сечении
km,
Затем вводится безразмерное время разработки.
τ – равное суммарному отбору жидкости в долях активных запасов.
(7)
где: q ж – дебит галереи по жидкости, м3/год
t – время, годы
QA – активные запасы.
Методом материального баланса в выражениях fH; 
(3,6,7) исключается переменная km и устанавливается зависимость между долей нефти fH и охватом 
и безразмерным временем разработки τ.
Эти зависимости и положены в основу дальнейших расчетов.
Для различных показателей неоднородности 
и соотношения вязкостей 
зависимости 
вычислены и табулированы в работе Баишева Б.П..
Начальные дебиты жидкости рядов скважин определяются по уравнениям интерференции Ю.П. Борисова, которые считаются постоянными на рассмотренном этапе и изменяются при переходе к следующему этапу, то есть после отключения ряда.
Система уравнений интерференции может быть представлена для полосовой залежи из трех рядов скважин (1нагнетательная и 2 добывающих) в следующем виде:
Где: 
– внешнее и внутрение сопротивление.
Линии тока
 | |||
 | |||
По зависимости 
или 
для всех значений τ находят величины
- доли нефти в потоке жидкости.
Затем определяются отборы нефти 
и воды 
– суммарный отбор нефти 
; воды 
; жидкости 
.
Кроме того метод Саттарова позволяет определить распределение скважин по дебитам и оценить их обводненость.
При этом принимается, что пласт состоит из определенного числа отдельных зон, а распределение проницаемости, по зонам подчиняется уравнению:
(8)
Где: ki – средняя проницаемость i-той зоны,
k п – параметрические распределения.
Этим же законом определяется и распределение проницаемости по мощности в каждой зоне.
Для установления темпа обводнения fi (τ i ) выбирается зона со средней проницаемостью k 0.
Для других зон безразмерное время определяется из соотношения:
Распределение скважин по дебитам находится по следующей формуле:
(9)
Где: q 0 – параметр распределения.
Затем по зависимостям f н (τ) или βн (τ) для всех значений τ выполняется расчет обводнения фонда скважин во времени.
Расчет ведется до тех пор, пока значение обводненности f в (τ) не достигнет заданного значения при отключении рядов скважин (98 – 99,9%)
ЛЕКЦИЯ
Расчет обводнения неоднородных пластов по методике института Гипровостокнефть (В.С. Ковалев, М.Л. Сургучев, Б.Ф. Сазонов)
Этот метод расчета процесса обводнения нефтяного пласта является дальнейшим развитием изложенных методов Ю.П. Борисова и М.М. Саттарова.
В методике института ГВН предусматривается более полный учет неоднородности коллектора (по проницаемости, пористости, нефтенасыщенности и коэффициенту вытеснения). Метод также предусматривает учет начальных водонефтяных зон, которые характеризуются параметром W.
Где: L 2 , L 1– расстояние от эксплуатируемой галереи (или ряда эксплуатирующих скважин) до внутреннего и внешнего контуров нефтеносности.
W изменяется от 0 до 1.
Если водо – нефтяная зона отсутствует (“запечатанные” залежи), то W =0 для залежи по всей площади нефтеносности подстилаемой водой (массивные или водоплавающие залежи) W =1. (0; 0,1; 0,2; 0,3; 0,4; 0,5; 0,6; ……1)
При расчете процесса заводнения нефтяного месторождения используется следующая схема:
 | |||
| |||
Пласт состоит из изолированных прослоев с различной проницаемостью, пористостью, начальной нефтенасыщенностью. При этом для каждого прослоя характерна своя величина коэффициента вытеснения нефти водой.
Истинная скорость перемещения фронта вытеснения, а следовательно и время обводнения каждого прослоя зависит от этих параметров.
где: 
mi – пористость,
Si – начальная нефтесыщенность,
η i– коэффициент вытеснения нефти водой.
При обводнении одинаковых по объему пропластков с различной эффективной пористостью из них можно получить не одинаковое количество нефти. При характеристике неоднородности продуктивных пластов рассматривается их комплексная неоднородность выраженная параметрами.
Где: ki – проницаемость i – го слоя.
mi , ρ i , η i – соответственно пористость, нефтенасыщенность и коэффициент
вытеснения i – го слоя.
Начальные продуктивные характеристики эксплуатации скважин находятся в результате их исследования при установившихся режимах фильтрации. Учитывается изменение коэффициентов фильтрации сопротивления A и B в результате деформации пласта коллектора, т.е. 
(коэффициенты увеличиваются), уравнение оттока воды для нагнетательной скважины имеет вид:
Где: q в.с. – расход воды в нагнетательной скважине.
(2)
Где: Рс – давление на забое нагнетательной скважины,
h в ; h г – обводненная и газонасыщенная толщина пласта,
μв, μг – соответственно вязкость воды и газа.
Из уравнения (2) находится расход воды q в. по одной “средней” нагнетательной скважине (при заданной депрессии ΔΡ в пласте).
Необходимое количество нагнетательных скважин равно:
(3)
Соответственно определяются другие показатели разработки 
на период ППД (поддержания пластового давления)
Зависимости между проницаемостью и другими параметрами пласта, установленные по изучению физико-геологических свойств продуктивных отложений Урало – Поволжья записываются в общем виде следующим образом:
k – проницаемость;
an , bn – постоянные коэффициенты, определяемые для каждого месторождения в результате исследования кернов по геофизическим данным.
Распределение параметра ω количественно оценивается коэффициентом вариации и описывается гамма – распределением, плотность которой имеет вид:
Где: υ;ωср. – параметры распределения (коэффициент вариации и среднее значение ω,
Г(в1+1) – гамма функция 
Для расчетов используется 15 распределений с коэффициентами вариации от 30,2 до 87,7%.
Расчеты также могут проводится и с применением логарифмически нормального закона распределения.
Плотность распределения которого имеет вид:
Где: σ; ωср – параметры распределения.
Соотношение для определения характеристик заводнения – доли нефти в потоке жидкости f н (τ) и β(τ) разработаны как для поршневого, так и с учетом непоршневого характера вытеснения.
Для поршневого вытеснения расчеты ведутся по характеристикам вытеснения с учетом различия вязкостей 
и скачкообразного изменения проницаемости в промытой зоне при условии наличия водо – нефтяной зоны.
Методика расчета полностью автоматизирована.
Расчеты проводятся в следующем порядке:
1. Изучается строение залежи и выбирается плотность сетки скважин в определенном диапазоне (например 400×400, 500×500, 600×600).
2. Определяется система разработки количество скважин, добывающих и нагнетательных для каждого из рассматриваемых вариантов.
3. Определяется средний дебит одной скважины по жидкости по результатам опробования скважин и гидродинамические исследования продуктивности залежи.
4. Рассчитывается параметр W.
5. Определяется соотношение вязкостей нефти и воды 
6. Выбирается закон распределения.
7. Рассчитываются активные запасы. 
8. По отдельной методике рассчитывается предельный дебит фонтанирования добывающих скважин.
9. Выбирается зависимость выбытия добывающих скважин при достижении предельной обводненности.
10. Принимается коэффициент компенсации объема добываемой жидкости закачкой. 
11. Расписывается ввод добывающих и нагнетательных скважин по годам.
12. Проводится расчет для каждого варианта по выбранной кривой до предельной обводненности продукции (99 – 99,9%).
13. Определяются основные технологические показатели и конечный коэффициент нефтеизвлечения.
Все рассмотренные методики расчета технологических показателей разработки нефтяных месторождений называются аналитическими.
Эти методики используются также для определения коэффициента нефтеизвлечения или нефтеотдачи. Прогнозирование роста обводненности продукции до 99,9% и определения соответствующих уровней добычи нефти, позволяет определить максимально возможное извлечение нефти, что составляет в итоге извлекаемые запасы залежи.
ЛЕКЦИЯ
Дата: 2018-12-28, просмотров: 276.
