СПОСОБ ОПРЕДЕЛЕНИЯ ЗАСТОЙНЫХ И СЛАБОДРЕНИРУЕМЫХ НЕФТЯНЫХ ЗОН В НИЗКОПРОНИЦАЕМЫХ КОЛЛЕКТОРАХ
DETERMINATION OF DEAD AND SLIGHTLY DRAINED OIL ZONES IN HETEROGENEOUS RESERVOIRS

Изобретение относится к нефтедобывающей промышленности, в частности к разработке нефтяных низкопроницаемых месторождений.
Известен способ определения местоположения застойных и слабодренируемых зон нефтяной залежи (пат. RU №2186204, опубл. 27.07.2002), осуществляемый путем выбора участка нефтяной залежи с добывающими и нагнетательными скважинами, определения для каждой добывающей скважины накопленного отбора нефти, жидкости и воды, времени эксплуатации скважины, построения карт изолиний, выделения на каждой из построенных карт слабодренируемых и застойных зон, сравнения выделенных зон на всех картах и принятия за слабодренируемые и застойные зоны нефтяной залежи участков в выделенных зонах, совпадающих на всех картах, отбором проб нефти из добывающих скважин на выбранном участке нефтяной залежи во времени не менее трех раз и не реже, чем через три месяца, определением по каждой пробе значения коэффициента светопоглощения нефти, построением карты изолиний на каждую дату отбора проб для значении коэффициента светопоглощения нефти и принятия участков, ограниченных изолиниями с максимальными, средними и минимальными значениями коэффициента светопоглощения нефти, соответственно, за активные, слабодренируемые и застойные нефтенасыщенные зоны нефтяной залежи.
Недостатком данного способа является его высокая трудоемкость, связанная с периодическими отборами проб нефти и определением коэффициента светопоглощения.
Наиболее близким техническим решением к предлагаемому изобретению является способ, основанный на выделении участков залежи, где образуется низкая скорость фильтрации пластовой нефти, обусловленная повышением вязкости нефти. Повышение вязкости нефти связывают с тем, что в процессе эксплуатации залежи происходит изменение свойств нефти при охлаждении и выделении из нее растворенного газа. Влияние вязкости на течение нефти наблюдается лишь при низких скоростях фильтрации, и в связи с этим выделяют критические градиенты давления - градиент динамического давления сдвига и градиент предельного разрушения структуры. При существовании в пласте фактических градиентов давления ниже градиента динамического давления сдвига и предельного разрушения происходит снижение охвата пласта фильтрацией, вследствие чего в пласте образуются слабодренируемые участки (В.В.Девликамов, З.А.Хабибуллин, М.М.Кабиров, «Аномальные нефти», Москва, «Недра», 1975 г., стр.135).
Для осуществления данного способа находят зависимость скорости фильтрации от градиента давления, из которой определяют критические точки: точка начала интенсивного роста объемной скорости фильтрации (характеризуется градиентом динамического давления сдвига) и точка перехода к линейному закону фильтрации (характеризуется градиентом давления предельного разрушения структуры). Градиенты давления определяют расчетным путем по средним значениям содержания структурообразующих компонентов нефти, величина которых не постоянна как по площади, так и по мощности. В связи с этим строят две карты распределения градиентов давлений: карту изоградиентов динамического давления сдвига и карту изоградиентов предельного разрушения структуры. Оценку расположения и размеров слабодренируемых зон производят сопоставлением значений характерных градиентов давлений с фактическими градиентами пластового давления, в связи с этим строят также карту фактических изоградиентов пластового давления. Границы слабодренируемых зон проводят путем совмещения карт расчетных характерных градиентов давления с картой распределения фактических градиентов давления. Участки, где фактические градиенты давления ниже градиента динамического давления сдвига и градиента давления предельного разрушения структуры, принимают за застойные и слабодренируемые.
Недостатком этого способа является то, что данный способ определения застойных и слабодренируемых зон основан на определении участков залежи, где образуются аномалии вязкости нефти, обусловленные изменением ее состава и свойств, проявляющиеся в основном в залежах, содержащих высокомолекулярные нефти и, соответственно, не приемлем для маловязких нефтей.
Основной решаемой задачей предлагаемого способа является применение способа для залежей маловязких нефтей (не обладающих аномалиями вязкости) в условиях низких проницаемостей. Технический результат от применения заявленного способа заключается в определении местоположения застойных нефтенасыщенных участков нефтяных низкопроницаемых залежей.
Поставленная задача решается тем, что на кернах проводят фильтрационные эксперименты со стационарным режимом фильтрации, строят зависимость скорости фильтрации жидкости от градиента давления, определяют значение предельного градиента давления, который соответствует изменению характера фильтрации жидкости. При этом дополнительно проводят на кернах различной проницаемости фильтрационные эксперименты с нестационарным режимом фильтрации, находят зависимость предельного градиента давления от проницаемости, строят карту модуля градиента пластового давления и карту проницаемости, наносят квадратную сетку на карту модуля градиента пластового давления и карту проницаемости, оценивают для каждой ячейки сетки значения модуля градиента пластового давления и коэффициента проницаемости, по зависимости предельного градиента давления от проницаемости вычисляют для каждой ячейки сетки значения предельного градиента давления, сравнивают полученные в ячейках значения предельного градиента давления со значениями модуля градиента пластового давления, выделяют ячейки, где модуль градиента пластового давления ниже значения предельного градиента давления, которые будут представлять собой застойную или слабо дренируемую зону нефтяной залежи. Способ осуществляется следующим образом.
1. Проводят серию лабораторных исследований по фильтрации жидкости через образцы кернового материала с различными проницаемостями. В каждом эксперименте проводят два типа фильтрационных экспериментов: стационарные и нестационарные. Нестационарная фильтрация используется для измерения скорости фильтрации при низких градиентах давления.
2. Строят график зависимости скорости фильтрации жидкости от градиента давления. На графике находят критическую точку, в которой происходит изменение угла наклона линии при переходе от высоких значений градиента давления к низким. Эта точка характеризует изменение режима фильтрации жидкости, а соответствующий ей градиент давления называется предельным градиентом давления. Нелинейный режим фильтрации проявляется при градиентах давления, меньших предельного градиента давления.
3. Строят график зависимости предельного градиента давления от проницаемости.
4. Строят карты проницаемости и пластового давления. На основе карты проницаемости и по найденной зависимости предельного градиента давления от проницаемости строят карту предельного градиента давления.
5. На основе карты пластового давления строят карту модуля градиента пластового давления.
6. На карты модуля градиента пластового давления и предельного градиента давления наносят квадратную сетку. Размеры сетки будут зависеть от необходимой степени точности получаемых результатов. Чем мельче сетка, тем более информативны результаты расчета на единицу площади залежи.
7. Внутри ячеек сетки определяют значения модуля градиента пластового давления и предельного градиента давления.
8. Сравнивают значения в соответствующих ячейках карты модуля градиента пластового давления и предельного градиента давления. Ячейки, в которых значения модуля градиента пластового давления меньше предельного градиента давления, будут являться застойными или слабодренируемыми участками залежи.
Пример.
2
3
3
Нефтяная залежь находится на стадии активной разработки. Выбираем краевой разбуриваемый участок залежи, включающий 153 добывающих и 78 нагнетательных наклонно-направленных скважин. На выбранном участке залежь характеризуется высокой расчлененностью и неоднородностью со следующими геолого-геофизическими параметрами: глубина залегания - 2600 м, эффективная нефтенасыщенная толщина - 28 м, средний коэффициент проницаемости по нефти - 0,0025 мкм, коэффициент пористости - 0,16, коэффициент нефтенасыщенности - 0,52, начальное пластовое давление - 25 МПа, вязкость нефти в пластовых условиях - 1,6 сП, плотность нефти в пластовых условиях - 870 кг/м, давление насыщения газом - 11,6 Мпа, газовый фактор - 70 м/т.
Для оценки застойных и слабодренируемых зон рассматриваемой разрабатываемой залежи необходимо получить зависимость предельного градиента давления от проницаемости. Для построения этой зависимости необходимо провести серию лабораторных исследований. В связи с чем, на керновых материалах, которые были отобраны из скважин, расположенных на рассматриваемом участке залежи, проводят лабораторные исследования по фильтрации жидкости. Выполняют 30 экспериментов, отличающихся проницаемостью керновых материалов. В свою очередь, в каждом эксперименте проводим два типа фильтрации.
0
0
Первый тип - стационарная фильтрация, когда при различных градиентах давления на керне измеряют скорость фильтрации флюида. Фильтрация в низкопроницаемом коллекторе не подчиняется линейному закону фильтрации Дарси. На фиг.1 показана характерная экспериментальная зависимость скорости фильтрации от градиента давления при проведении экспериментов по стационарной фильтрации. На фиг.1 можно выделить линейный участок, причем его экстраполяция отсекает на оси абсцисс эффективный начальный градиент давления (Н) и участок в области низких градиентов давления. В результате обработки экспериментов по стационарной фильтрации можно определить угол наклона линейной части в зависимости скорости фильтрации от градиента давления, тангенс угла которого характеризует проницаемость линейной части, и эффективный начальный градиент (Н).
Эксперименты по стационарной фильтрации не позволяют измерить скорость при достаточно низких градиентах давления. Поэтому проводят второй тип экспериментов - нестационарная фильтрация. В данном случае после установления стационарного состояния перекрывают вход и выход на образце фильтруемого керна и измеряют изменение градиента давления во времени. Для обработки экспериментов по нестационарной фильтрации используют численное решение уравнения пьезопроводности на давление с соответствующими граничными условиями:
<mrow><mfrac><mo>∂</mo><mrow><mo>∂</mo><mi>x</mi></mrow></mfrac><mo stretchy="false">(</mo><mi>χ</mi><mfrac><mrow><mo>∂</mo><mi>P</mi></mrow><mrow><mo>∂</mo><mi>x</mi></mrow></mfrac><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><mo>∂</mo><mi>P</mi></mrow><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac></mrow>
<mrow><mfrac><mrow><mo>∂</mo><mi>P</mi></mrow><mrow><mo>∂</mo><mi>x</mi></mrow></mfrac></mrow>
(x=0; f>0)=0
<mrow><mfrac><mrow><mo>∂</mo><mi>P</mi></mrow><mrow><mo>∂</mo><mi>x</mi></mrow></mfrac></mrow>
(x=L; t>0)=0
<mrow><mtext>P(x; t</mtext><mo>=</mo><mtext>0)</mtext><mo>=</mo><msub><mrow><mtext>(P</mtext></mrow><mtext>2</mtext></msub><msub><mrow><mtext>-P</mtext></mrow><mtext>l</mtext></msub><mtext>)*</mtext><mfrac><mtext>x</mtext><mtext>L</mtext></mfrac><mo>+</mo><msub><mi>P</mi><mn>1</mn></msub></mrow>
где коэффициент пьезопроводности:
<mrow><mi>χ</mi><mo>=</mo><mfrac><mrow><mi>k</mi><mo>*</mo><msub><mi>k</mi><mi>r</mi></msub></mrow><mrow><mi>μ</mi><msub><mi>c</mi><mi>t</mi></msub></mrow></mfrac></mrow>
,
r
t
1
2
k - абсолютная проницаемость, k - относительная фазовая проницаемость, µ - вязкость, c - общая сжимаемость системы, L - длина образца керна, P и P - левое и правое значение давления на концах керна в нулевой момент времени.
При интерпретации экспериментов по нестационарной фильтрации решается обратная задача восстановления проницаемости в области низких градиентов давления. В качестве невязки используют ошибку между решением математической модели и результатом эксперимента. На фиг.2 показан пример результатов экспериментов нестационарной фильтрации (светлая линия с маркерами) - зависимость градиента давления от времени. Здесь можно четко выделить две зоны: быстрое падение градиента в начальный момент (Зона 1) и последующее медленное изменение динамики падения (Зона 2). Первая зона соответствует линейному участку на фиг.1, а вторая зона - области низких градиентов давления. Сплошной темной линией без маркеров на фиг.2 показаны результаты численного решения обратной задачи по восстановлению проницаемости в области низких градиентов давления (на фиг.2 совпадает с кривой результатов эксперимента нестационарной фильтрации), которые позволяют в области низких градиентов давления определить угол наклона зависимости скорости фильтрации от градиента давления, тангенс угла которого характеризует проницаемость в области низких градиентов давления.
0
L
В результате совместного проведения экспериментов по стационарной и нестационарной фильтрации можно определить проницаемость линейной части, эффективный начальный градиент (Н) и проницаемость в области низких градиентов давления. Используя данные параметры, можно определить предельный градиент давления (H). На фиг.3 показан предельный градиент давления, который получается путем пересечения двух линий: первая в области низких градиентов давления, а вторая в области высоких градиентов давления с эффективным начальным градиентом давления. Темная сплошная линия характеризует зависимость скорости фильтрации от градиента давления при стационарной фильтрации, а пунктирная светлая линия характеризует численную интерпретацию результатов эксперимента по нестационарной фильтрации. Предельный градиент давления характеризует изменение характера фильтрации.
L
L
Для каждого образца кернового материала, при каждом опыте получают зависимость, подобную зависимости на фиг.3, с точкой изменения характера фильтрации - предельным градиентом давления (H). По результатам выполнения всех 30 экспериментов получают зависимость предельного градиента давления (H) от абсолютной проницаемости (фиг.4).
L
L
L
L
Далее, по результатам расчетов на гидродинамической модели, имеющей шаг сетки 50 м, строят карту текущего пластового давления. На фиг.5 представлена карта текущего пластового давления рассматриваемого участка залежи, где цвет в ячейках обозначает значение текущего пластового давления, величину которого можно узнать по соответствующему цвету на шкале, расположенной справа от карты текущего пластового давления. Протяженные осветленные участки на карте - это участки высокого давления в окружении трещин авто-ГРП, которые образуются на нагнетательных скважинах при забойных давлениях, превышающих давление разрыва пласта. На основе карты текущего пластового давления строят карту модуля градиента пластового давления (фиг.6). На фиг.6 темные участки представляют низкие градиенты пластового давления. Также для рассматриваемого участка строят карту проницаемости (фиг.7). На основе карты проницаемости и по найденной зависимости предельного градиента давления H от проницаемости (фиг.4) строят карту предельного градиента давления H (фиг.8). Сравниваем карты модуля градиента пластового давления (фиг.6) и предельного градиента давления H (фиг.8). Ячейки, в которых значения модуля градиента пластового давления меньше предельного градиента давления H, будут являться застойными участками залежи. На фиг.9 представлена карта размещения застойных или слабодренируемых зон, где светлые участки являются застойными или слабодренируемыми зонами, а темные - незастойными зонами. Таким образом, предлагаемый способ позволяет определить местоположение застойных или слабодренируемых нефтенасыщенных участков низкопроницаемых залежей вне зависимости от вязкости нефтей.